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1.  INTRODUCTION 

This  document  is  a  preliminary  report  on  a  survey  of  the 
factors  affecting  the  performance  of  a  filter.  Although  specifi¬ 
cally  directed  at  optical  filters,  most  of  the  discussion  is 
quite  general,  applying  to  the  propagation  of  any  wave  in  a 
medium  with  a  spatially  modulated  response  or  susceptance.  The 
intent  of  this  work  is  to  examine  the  factors  affecting  filter 
performance  with  a  view  toward  perfecting  filter  design.  This 
effort  began  with  a  survey  of  the  optical  properties  of  materials 
which  could  be  used  in  a  filter.  The  assumption  was  that  ulti¬ 
mate  filter  performance  is  primarily  limited  by  materials 
questions.  During  this  survey  the  author  quickly  became  aware  of 
the  need  to  develop  a  rating  system  for  material  properties, 
based  on  the  relationships  between  filter  performance  and  the 
susceptance  of  the  various  parts  of  the  filter.  Thus  it  became 
necessary  to  first  establish  this  relationship,  a  task  which 
forms  the  main  body  of  this  work. 

1 . I  Design  Needs 

The  filter  design  engineer  is  severely  hampered  if  he  does 
not  have  the  freedom  to  specify  the  susceptance  (this  implies  the 
refractive  index)  of  the  medium  at  a  given  location.  A 
dielectric  can  be  tailored  to  an  arbitrary  specification  between 
the  limiting  value  of  the  dielectric  constants  of  a  given  set  of 
materials  by  making  an  appropriate  solution  of  these  materials. 
Unfortunately  there  is  no  ready  formula  for  the  ideal  solution 
(this  question  is  addressed  in  Paragraph  2.1)  and  no  solution 
will  be  perfectly  homogeneous.  The  impact  of  this  fact  on  wave 
propagation  is  discussed  in  Paragraph  2.2. 
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Bragg  filters  are  based  on  the  interference  effects  upon 
wave  propagation  due  to  spatial  variations  in  the  susceptance  of 
the  medium.  Traditional  optical  filters  have  modulated  the 
refractive  index  profile  in  discrete,  discontinuous  steps.  The 
index  response  between  these  steps  is  flat.  This  method  achieves 
great  simplicity  in  both  manufacture  and  design.  The  motivation 
for  this  work  is  the  production  of  a  new  type  of  filter--the 
graded  index  filter.  We  define  a  graded  index  filter  as  one  in 
which  the  medium  susceptance  (and  therefore  the  dielectric 
constant)  is  smoothly  and  continuously  varying,  i.e.,  one  in 
which  there  are  no  sharp  discontinuities  or  flat  regions  in  the 
susceptance.  Currently  a  sinusoidal  modulation  is  popular  in 
des  ign. 

Much  of  the  impetus  for  the  development  of  graded  index 
filters  is  based  upon  the  speculation  that  the  freedom  to  con¬ 
tour  the  index  profile  at  will  should  allow  the  design  engineer 
extra  flexibility  in  the  design  of  filters  with  ultra  narrow 
stop,  or  reflectance  bands,  the  location  of  several  such  bands  at 
arbitrary  frequencies,  and  the  avoidance  of  absorption  and/or 
scattering  associated  with  material  discontinuities.  The  author 
has  uncovered  nothing  in  this  study  to  support  these  beliefs.  On 
the  contrary,  the  analysis  presented  here  suggests  that  details 
of  the  modulation  profile  have  little  bearing  on  filter  perfor¬ 
mance  in  the  frequency  range  for  which  most  simple  filters  are 
designed . 

Wave  propagation  in  a  periodic  medium  is  governed  primarily 
by  fundamental  considerations  without  regard  to  the  details  of 
the  system  in  question  (some  of  these  fundamental  features  are 
discussed  in  the  Appendix).  Consequently,  a  wide  range  of 
systems  which  are  usually  treated  as  unique  share  important 
features  with  optical  filters.  Examples  are  acoustic  filters 
(mufflers),  electronic  energy  bands  in  solids,  lattice  dynamics, 


transmission  line  theory,  and  waveguide  theory.  Insight  to  the 
optical  filter  problem  can  be  gained  by  reviewing  the  large  body 
of  work  developed  in  these  other  related  areas.  In  particular, 
it  is  universally  found  that  a  periodic  modulation  of  the  medium 
response  will  produce  not  just  a  single  stop,  or  reflectance 
band,  but  a  series  of  such  bands  with  increasing  frequency.  In 
solid-state  band  theory,  these  stop  bands  are  referred  to  as  band 
gaps.  The  higher  order  band  structure  is  related  to  the  details 
of  the  modulation  profile.  In  the  event  a  filter  is  designed  for 
a  multi-octave  frequency  range,  the  details  of  the  index  contour 
will  become  important. 

1 . 2  Approach 

The  analysis  adopted  in  the  following  review  of  filter  per¬ 
formance  is  analytic.  The  philosophy  is  that  an  analytic  (in 
contrast  to  a  numerical)  approach  has  advantages  for  conferring 
insight  on  the  essential  physical  parameters  that  determine 
filter  performance.  In  the  process  of  pursuing  this  approach  a 
second  goal  developed:  the  application  of  analysis  to  optical 
filter  design,  especially  the  design  of  continuous  gradient 
filters.  A  general  procedure  for  designing  a  filter  is  given  in 
Figure  1.1,  which  is  taken  from  a  reference  on  the  design  of 
filters  in  transmission  lines.  1  This  report  begins  to  lay  the 
groundwork  for  development  of  an  optimum  transfer  function  and  a 
determination  of  appropriate  values  for  the  filter  parameters. 

The  author  believes  that  this  analysis  will  answer  questions  on 
the  effect  of  discontinuities  and  of  continuously  varying  suscep¬ 
tibilities  on  filter  performance,  and  the  degree  to  which  an 
ideal  notch  filter  can  be  approximated. 

An  analytic  approach  has  special  advantages  when  a  problem 
becomes  complex.  An  example  is  the  extension  of  simple,  one¬ 
dimensional  filters  with  purely  real  refraction  indices  to  three- 


THE  FLOWCHART  OF  FILTER  DESIGN  OPERATIONS 


Figure  1. 


Approach  to  filter  design  used  in  transmission  line 
theory  ^ 
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dimensional  filters  with  complex  susceptibilities.  This  problem 
can  be  handled  analytically  by  a  resolvant,  or  Green's  function 
approach.  The  Green's  function  method  is  a  generalization  of  the 
concepts  addressed  in  this  report. 

Section  3  presents  an  overview  of  filters,  of  different 
means  of  describing  their  performance,  and  of  the  diverse  sub¬ 
jects  which  share  essential  features.  The  critical  conclusion  of 
Section  3  is  that  a  complete  description  of  the  performance  of  a 
filter  may  be  obtained  once  the  wave  equation  is  solved  over  the 
fundamental  repeat  region  (unit  cell)  of  the  filter  cascade. 

The  solution  of  the  wave  equation  is  then  discussed  in  Section  4, 
with  special  emphasis  on  the  two  cases  of  special  interest,  step 
function  and  sinusoidal  modulation.  These  two  cases  are  com¬ 
pared,  and  some  general  observations  are  made  about  the  behavior 
of  the  solutions,  particularly  the  appearance  of  stop  and  pass 
bands . 

Sections  5  and  6  discuss  various  means  by  which  the  solution 
of  the  wave  equation  may  be  converted  to  an  algebraic  problem. 
Algebraic  procedures  are  especially  useful  in  describing  multiple 
periodic  filters,  as  well  as  terminations,  junctions,  and  other 
local  departures  from  the  periodic  profile.  They  are  a  prime 
analytic  tool.  The  mathematical  basis  for  Section  5  is  given  in 
Appendix  A. 

Several  candidates  for  optical  filter  materials  are 
discussed  in  Section  7,  largely  in  terms  of  crystal  structure  and 
lattice  parameters,  since  these  must  be  well  matched  to  obtain  a 
minimal-defect  solution.  We  call  attention  to  the  as-yet-little- 
studied  layered  materials,  in  particular  the  transition  metal 
chalcogenides.  Our  conclusions  and  suggestions  for  future  work 
are  given  in  Section  8.  Although  this  report  is  largely  the 
result  of  an  extensive  literature  review,  every  section  contains 
significant  original  contributions  by  the  author. 


2. 


PRELIMINARY  CONSIDERATIONS  ON  DIELECTRIC  BEHAVIOR 


The  following  three  sections  will  give  a  mathematical 
description  of  light  propagation  in  an  ideal  medium.  In  this 
section  we  comment  on  formulas  appearing  in  the  literature  for 
tailor-made  dielectrics  and  on  the  effect  of  a  material's  imper¬ 
fections  on  its  ideal  behavior. 

2.1  Dielectric  Constant  of  a  Mixture 


Before  we  examine  how  departures  from  an  ideal  profile  will 
affect  the  response  of  an  inhomogeneous  medium,  we  should  first 
obtain  some  insight  into  the  dielectric  behavior  of  a  solution. 

To  our  knowledge,  this  matter  has  not  been  adequately  explored. 

The  dielectric  constant  of  a  material  is  a  macroscopic  property 
with  a  nonlinear  dependence  on  the  molecular  (microscopic) 
polarizability.  The  reason  for  this  nonlinearity  is  the 
depolarizing  response  to  an  external  field  created  at  any  discon¬ 
tinuity  in  the  susceptibility  of  the  medium.  Consequently, 
accurate  models  of  dielectric  behavior  must  take  account  of  the  effect  of 
the  surrounding  medium  on  the  microscopic  susceptibility  of  a 
material.  Unfortunately,  this  behavior  is  very  complex.  We  now 
derive  a  general  result  for  a  multi-component  dielectric  mixture. 

The  relation  of  our  result  to  earlier  work  and  suggestions  for 
extensions  follow. 

We  denote  the  applied  electric  field,  the  local  effective 
electric  field,  and  the  net  polarization  field  respectively  by  E, 

E',  and  P.  We  will  model  the  medium  by  a  solution  of  molecules 
with  ellipsoidal  shapes.  The  polarizability  of  a  molecule  of  a 
given  species  will  be  denoted 


We  have 


(2.2) 


or 
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(2.10 


These  relations  neglect  permanent  dipole  moments.  In  fact, 
most  materials  of  interest  in  hard  filters  are  polar;  neverthe¬ 
less  our  neglect  of  permanent  dipole  moments  remains  valid  at 
optical  frequencies.  In  the  intermediate  infrared  region  the 
permanent  molecular  dipoles  can  respond  to  the  applied  field,  so 
we  must  make  the  substitution 


»i  +  /(3kT  +  b  n.  ) 


(2.11) 


for  dipoles  of  series  i  with  concentration  ni  at  a  tem¬ 
perature  T.  In  Equation  (2.11),  b  is  an  empirical  constant. 

We  obtain  a  relation  between  the  macroscopic  parameters  e 
and  £i  and  t2  for  a  two-component  solution  of  molecular  densities 
nji  and  n2  by  an  iteration  of  Equation  (2.8): 
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This  yields 
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which  leads  to 


”l^l^l  ”2^2^2 

"1^1  ^  "2^2 


(2.14) 


where 


l/(Ae.  +  1  -  A} 


(2.15) 


a 


When  A  =  1/3  and  0  we  respectively  obtain  the  relations  for  the 
Clausius-Mossotti  and  "Drude,"  or  more  properly,  Sellmeier  models 
given  by  Bottcher^  and  Jacobsson. 2  We  note  that  the  simplified 
form  found  in  Equation  (2.14)  is  only  valid  if  A  is  the  same  for 
both  species.  It  can  be  extended  to  include  a  summation  over 
several  species  in  a  solution. 


We  can  extend  Equation  (2.14)  by  induction,  and  give  here  an 
illustration  of  the  first  step  in  the  induction  argument,  a  three- 
component  solution.  We  consider  the  solution  to  be  composed  of  a 
two-component  solution  plus  the  third  component,  and  use  the 
subscript  o  to  refer  to  the  two-component  solution.  Then  from 
Equation  (2.14)  we  have 


e  = 


n  a  c  +  n-a.e^ 
OOP _ 333 

n  +•  n-a, 
o  3  3 


(2.16) 


Again  from  Equation  (2.14),  we  rewrite  Equation  (2.16)  as 


(n^+n2J  ^1  ^"^”2  ^  ^  ^  3®  3^  3 


e  = 


'^1®1^"2®2 


"^3^3 


(2.17) 


Now 


^o  =  ^^"^o  +  1  -  A) 


-1 


(2.18) 


or 


a  =  1/[A  -A-LJ  + 

o  nj^aj^+n2a2 


1  -  A' 


(2.19) 


Consequently 


"l  "2 

(ni  +  n2)a^  - 


(2.20) 
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From  Equation  (2.15),  we  see  that 

(Hj^  +  n2)  =  n^aj^  +  n2a2  .  (2.21) 

Substitution  of  Equation  (2.21)  in  Equation  (2.17)  gives 

3 

I  n.a.e. 

e  =  — ^ -  .  Q.E.D.  (2.22) 

I  n.a. 

Clearly,  we  may  iterate  this  argument  to  include  any  number  of 
species  of  the  same  A  coefficient. 

In  general,  e,  x»  and  a  will  be  complex.  Equation  (2.14)  is 
valid  for  complex  dielectric  constants,  as  may  be  seen  from  its 
derivation.  Part  of  its  utility  is  that  it  gives  a  result  for 
complex  dielectric  constants  without  reference  to  the  details  of 
the  microscopic  polarizability.  From  our  derivation,  we  see  that 
Equation  (2.14)  should  also  include  a  reasonable  amount  of 
dispersion,  by  letting  e  -*■  e(a)). 

In  Equation  (2.14)  we  recover  the  value  for  a  pure  material 
if  either  nj^  or  n2  0,  as  we  should.  In  a  mixture  t  is  reduced 
from  a  value  given  by  a  linear  interpolation  between  and  e2 
for  A  ^  0  due  to  the  interaction  between  molecules.  To  see  this, 
consider  a  solution  when  nj  =  n2.  Then  we  find 

2Ae  £2  +  (1  -  A)  [e,  +  E  ) 

e  =  - — - - - —  .  (2.23) 

2  ( 1  -  A]  +  a(  Ej^  +  E2  ) 

Comparing  this  with  (ej^  +  t2]/2  we  obtain  the  variation  in  e 
from  its  value  found  from  a  linear  interpolation,  6e: 

A( E,  -  E2 

<Se  =  -  .  (2.24) 

4[2(1  -  A)  +  A(e^  +  £2)] 
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We  see  that  6e  will  increase  with  both  A  and  £2  -  £i.  This  is 
consistent  with  experiment  and  the  discussion  given  by 
Jacobsson. 2  However,  part  of  Jacobsson's  treatment  contains 
errors,  since  it  is  based  on  a  faulty  premise  which  we  now 
descr  i  be . 

If  an  ellipsoid  with  dielectric  constant  £2  is  embedded  in  a 
uniform  medium  of  dielectric  constant  ej  then  a  net  simple  dipole 
field  proportional  to  £2  -  will  be  added  to  the  effective 
field  that  would  occur  in  the  absence  of  the  ellipsoid.^  Using 
this  result.  Polder  and  Van  Santen^  obtained  the  effective 
dielectric  constant  for  a  dilute  suspension  of  such  spheres. 

£  =  £1  +  Cei(£2  -  ei)/[Ae2  +  (l  "  A)ei]  ,  (2.25) 

where  C  is  the  volume  concentration  of  the  ellipsoids.  This 
result  is  readily  extended  to  include  ellipsoids  of  varying 
dielectric  constant  and  shape.  A  similar  result  is  found  by 
Landau  and  Lifshitz  for  spheres. ^  However,  these  authors  were 
mindful  of  the  difference  between  a  dilute  suspension  and  a 
strong  solution,  noting  that  Equation  (2,25)  can  only  be  taken  as 
an  approximation  for  a  solution  at  low  concentrations.  We  may 
attempt  to  extend  Equation  (2.25)  by  using  the  approximation 

£  =  £1  +  CeNe2  -  ei)/[Ae2  +  (l  -  A)e]  ,  (2.26) 

as  Polder  and  Van  Santen  in  fact  did,  but  this  formula  also 
breaks  down  as  C  increases,  as  may  be  seen  by  comparison  with 
Equation  (2.12)  using  C  =  x\2/(n\  +  n2).  Jacobsson  ignores  this 
limitation. 

Equation  (2.26)  is  a  macroscopic  result.  It  pertains  to 
particles  in  a  homogeneous  background.  We  may  well  treat  the 
interaction  between  a  molecule  and  distant  molecules  by  an 
average  field, ^  but  must  be  careful  about  our  treatment  of  local 
interactions.  Jackson, ^  shows  that  the  Clausius-Mossotti  result 
(A  =  1/3)  should  work  well  for  distant  interactions,  with  local 
fields  accounted  for  by  a  perturbation  parameter  s. 


A  =  1/3  +  s  . 


(2.27) 


For  synunetrical  systems,  such  as  cubic  crystals  and  pure 
amorphous  materials,  the  various  contributions  which  include  s 
cancel.  This  is  why  the  Claus ius-Mossott i  equation  is  so  useful. 

Departures  observed  from  the  values  predicted  using 
A  =  1/3  are  observed  in  bicomponent  films. ^  This  may  be  due  to 
a  lack  of  local  isotropy,  but  could  well  be  due  to  local  depar¬ 
tures  from  stoichiometry  and  voids.  It  is  impossible  to  say 
without  a  thorough  characterization  of  the  films.  As  can  be  seen 
from  Equation  (2.25),  voids  can  have  a  dramatic  impact  on  the 
dielectric  constant. 

The  derivation  leading  to  Equation  (2.8)  is  inexact.  A  more 
accurate  formu'  tion  can  be  obtained  following  the  Onsager  devel¬ 
opment  for  sph»  es.l  In  fact,  we  could  use  a  statistical  treat¬ 
ment.  1  In  our  treatment  we  have  swept  statistical  variations  in 
the  local  dielectric  environment  of  a  molecule  under  the  rug, 
burying  the  resulting  parameter  s  in  our  parameter  A.  In  other 
words,  our  treatment  is  intended  only  as  a  model  that  will 
include  the  effect  of  an  inhomogeneous  local  environment. 

However,  we  believe  it  gives  a  better  picture  than  an  uncritical 
use  of  Equation  (2.26).  It  cannot  readily  be  extended  to  include 
solutions  of  materials  with  different  values  for  A.  However, 
Equation  (2.8)  should  be  adequate  for  our  purposes.  In  fact, 
until  it  is  established  that  departures  from  the  Claus ius-Mossott i 
relation  are  in  fact  due  to  molecular  shape  factors  rather  than 
film  preparation,  we  may  well  choose  A  =  1/3.  For  design  pur¬ 
poses  it  is  probably  best  to  use  empirical  values  for  the 
dielectric  constants  of  mixtures,  since  they  are  so  heavily 
dependent  on  fabrication  processes. 
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2 . 2  Rayleigh  Scattering  from  I nhomoqene i t ies  in  e 


Local  variations  from  the  desired  dielectric  coefficient 
will  mar  its  performance.  We  are  concerned  here  with  random, 
local  imperfections,  rather  than  regular  one-dimension  discon¬ 
tinuities  in  the  dielectric  profile.  These  imperfections  will 
include  voids,  departures  from  stoichiometry  due  to  lack  of 
control  during  fabrication,  departures  from  stoichiometry  due  to 
compositional  segregation  following  fabrication,  surface  rough¬ 
ness  at  film  interfaces,  impurities,  density  fluctuations  in 
amorphous  films,  the  occurrence  of  microcrystals  in  an  amorphous 
matrix,  and  polycrystal 1 ine  films  containing  more  than  one  phase, 
or  randomly  oriented  crystals  with  anisotropic  refractive  indi¬ 
ces.  We  assume  that  all  these  imperfections,  if  they  occur,  will 
be  small  in  size.  This  will  allow  us  to  make  a  Rayleigh  scat¬ 
tering  approximation.  Consequently  we  will  not  need  to  consider 
the  shape  of  the  scattering  centers,  although  except  for  the  case 
of  voids  we  may  easily  extend  our  discussion  to  include  centers 
of  any  shape. ^  Similarly,  we  will  ignore  gradations  in  the  index 
profile  of  imperfections,  and  model  them  by  spheres  of  uniform 
dielectric  strength.  A  more  accurate  description  will  add  little 
to  our  conclusions  at  the  cost  of  great  complications,  since  we 
can  account  for  most  effects  by  choosing  an  effective  radius,  a, 
for  the  scattering  centers.  We  will  also  neglect  coherence 
effects.  We  denote  the  scattering  cross  section  by  Og  and  the 
absorption  cross  section  by  03.  We  let  the  local  dielectric 
constant  of  the  medium  be  em?  and  that  of  the  imperfection  be  cj. 
We  have 


and 


~  T-  ^ 


a  =  4irk 

d 


(2.28) 


(2.29) 


(2.30) 


assuming  that  the  imaginary  part  of  e,  e"  is  small  compared  to 
its  real  part,  e'.  Here 

k  =  2Tr/e^/XQ,  (2.31) 

is  the  wave  vector  in  the  medium  surrounding  the  scattering 
center  corresponding  to  a  freespace  wavelength  Xq*  We  see  that 


e  .•  +  2  I 
1  m  J 


(  ka ) 


II  11 


(2.32) 


Since  we  are  considering  a  regime  in  which  ka  <<  1,  absorp¬ 
tion  will  be  important,  and  in  fact  will  probably  dominate.  An 
exception  to  this  rule  will  again  be  the  case  of  voids. 


From  Equation  (2.28)  we  see  that  scattering  will  be  far  more 
severe  in  the  optical  region  than  in  the  infrared.  Absorption 
will  also  be  more  severe,  although  not  proportionately.  We  will 
assume  that  the  density  of  scattering  centers  is  sufficiently 
low,  so  that  we  may  restrict  ourselves  to  single  particle  scat¬ 
tering.  The  intensity  at  a  depth  z  in  the  filter,  I(z),  will  be 
given  by 


I(z)  =  Iq  exp[/dz  n(z)  Oe(z),  (2.33) 

where  Iq  is  the  input  intensity  and 

<30—  .  (2.34) 

We  see  from  Equation  (2.29)  that  the  total  absorption  per  unit 
length  due  to  a  given  type  of  imperfection  is  proportional  to  the 
fractional  volume  of  scattering  centers  in  the  material.  Thus  for 
small  centers  the  relative  importance  of  a  given  type  of  imperfection 

II  M 

IS  roughly  proportional  to  C  (?[  -  e^,)  r  where  C  is  the  fractional 
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volume  of  that  type  of  impurity 


We  now  present  models  for  the  dielectric  behavior  of  several 
imperfections.  Voids  will  have  the  strongest  scattering  for  a 
given  size  scattering  center,  but  will  not  contribute  to  absorp¬ 
tion.  We  have 


1  , 


(2.35) 


A  tendency  towards  separation  or 
will  drive  both  and  toward 
complete  separation  we  will  have 


decomposition  of  the  solution 
ei  or  £2*  the  limit  of 


ei  £2  *  Em  ei  .  (2.36) 

We  may  treat  a  partial  decomposition  in  the  same  manner  as  a 
variation  in  stoichiometry.  Let 

e  =  Cl  +  n2 £2  .  (2.37) 

Then  for  a  fractional  change  6  in  composition  we  will  have 

e  =  (niei  +  n2e2)  +  '^{ei  -  e2  )  ,  (2.38) 


so  that 


ei  =  6(ei  -  €2) 


(2.39) 


and 


*('1  -  tj)  •  (2.40 

We  expect  6  to  be  on  the  order  of  a  few  percent. 

If  the  dielectric  constant  changes  due  to  a  change  in  phase, 
or  in  the  case  of  an  anisotropic  medium,  to  a  change  in  orien¬ 
tation,  we  model  it  as 
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(  2.41  ) 


so  that 


"i  -  "m  =  1  •  ^2.42) 

We  expect  that  A  may  be  on  the  order  of  several  percent.  There 
may  also  be  a  small  change  in  the  dielectric  constant  at  the 
boundaries  of  grains  in  polycrystal ine  material  or  at  film  inter¬ 
faces.  Here  also  we  expect  a  relations  like  those  of  Equations 
(2.41)  and  (2.42),  with  the  fractional  volume  of  scattering  centers 
proportional  to  the  surface  to  volume  ratio  of  the  grains  or 
f i 1ms . 
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3. 


WAVE  PROPAGATION  IN  A  PERIODIC  STRUCTURE 


A] I  periodic  structures  act  as  filters  to  waves  passing 
through  them,  i.e.,  they  have  pass  and  stop  bands.  The  problem 
we  wish  to  solve,  the  optical  filtering  properties  of  a  material 
with  a  modulated  index  of  refraction,  is  very  similar  to  other 
problems  in  acoustics , electrical  transmission  lines, 6, 7 
microwave  wavegu  ides ,  solid-state  band  theory,  lat¬ 
tice  dynamics .  We  can  gain  a  great  deal  of  insight  from  the 

work  done  in  these  related  areas.  The  first  thing  that  becomes 
apparent  in  comparisons  with  these  other  systems  is  that  the  stop 
and  pass  bands  do  not  depend  greatly  on  the  specific  nature  of 
the  wave  in  question;  rather  they  simply  arise  from  the  regular 
variation  of  the  response  of  the  transmitting  medium  to  the 
disturbance  of  the  wave. 

Procedures  for  analytic  treatment  of  problems  of  this  kind 
are  well  developed.  We  now  summarize  them  and  give  examples.  We 
will  then  comment  on  the  possibilities  for  inversion,  investi¬ 
gating  high-order  band  behavior,  two  and  three  dimensional  calcu¬ 
lations,  and  the  similarities  and  differences  between  the 
response  of  filters  with  slightly  different  index  profiles. 

3 . 1  Graphic  Displays  of  Performance 

The  result  of  a  calculation  will  be  a  curve.  It  may  be  a 
transmission  or  reflection  coefficient  curve, 44,45  displaying  the 
transmission  coefficient,  T,  or  reflection  coefficient,  R,  versus 
the  incident  wavenumber,  ko,  it  may  also  be  a  dispersion 
curve, 46-50  displaying  the  real  part  of  the  frequency  of  the 
incident  wave,  oi,  (which  is  equivalent  to  kQ)  versus  the  wave- 
number  in  the  medium,  k.  These  two  types  of  curves  contain  simi¬ 
lar  information,  since  attenuated  transmission  is  accounted  for 
by  a  complex  value  for  u,  and  we  can  easily  plot  the  imaginary 
part  of  u)  as  well  as  its  real  part.  When  the  attenuation  is 
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small,  there  are  some  advantages  to  using  the  dispersion  curves 
instead  of  reflection  curves.  The  dispersion  curves  give  a  clear 
picture  of  the  pass  and  stop  bands  above  the  fundamental  bands. 
They  give  a  good  picture  of  behavior  in  three-dimensional 
systems,  are  very  useful  in  analyzing  mode  coupling  and  the  den¬ 
sity  of  states,  or  channel  carrying  capacity  in  a  given  frequency 
range,  Au.  Consequently  dispersion  curves  are  commonly  used  in 
solid-state  physics  and  waveguide  transmission  problems. 
Apparently  this  type  of  analysis  has  been  relatively  little  used 
in  optics.  Optical  problems  have  more  commonly  been  expressed  in 
terms  of  reflection  and  transmission  coefficients,  which  have 
advantages  when  absorption  by  the  medium  is  large  or  when  the 
ratio  of  the  filter  thickness,  t,  to  the  length  characterizing 
its  periodicity  (unit  cell  length),  d,  is  small.  Since  disper¬ 
sion  curves  are  useful  in  obtaining  physical  insight,  we  will 
refer  to  dispersion  curves  in  most  of  our  discussion. 

There  are  two  basic  starting  points  for  solving  wave 
transmission  problems.  The  first  is  a  direct  solution  of  the 
partial  differential  equation  (PDE)  for  the  wave  amplitude  and 
phase,  i(/.  The  second  is  a  conversion  from  a  PDE  to  a  difference 
equation.  The  difference  equation  approach  is  most  useful  when 
the  index  modulation  is  large  over  short,  well  separated  regions 
(tight-binding  approximation)  or  when  the  modulation  is  small 
(nearly  free  case).  In  the  latter  case  the  difference  equation 
is  obtained  by  a  Fourier  transformation  of  the  PDE.  Of  course, 
the  direct  solution  of  the  PDE  is  most  useful  when  the  index 
modulation  over  a  given  interval  has  a  behavior  close  to  that  in 
which  solutions  are  known.  Perturbation  theory  may  be  used  with 
either  approach. 

Each  type  of  analysis  is  most  useful  when  expressions  for 
the  dispersion/attenuation  or  reflection/transmission  curves  are 
given  in  terms  of  the  functional  form  of  the  index  modulation. 
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However,  a  good  qualitative  understanding  of  a  problem  is  often 
possible  in  terms  of  integrals  which  may  be  obtained  numerically 
for  specific  instances.  Figure  3.1  shows  the  relationships  bet¬ 
ween  various  functions  commonly  used  in  different  approaches  to 
the  problem. 

The  solution  to  the  wave  equation  and  its  derivative  may  be 
thought  of  as  a  vector,  so  that  the  quantity  relating  the  input 
and  output  values  of  this  vector  is  a  matrix.  In  an  optics 
problem  we  actually  need  to  know  four  functions  to  completely 
specify  one  of  these  vectors,  (e.g.,  for  propagation  in  the  x 
direction  we  need  to  know  Ey »  and  Hy ) ;  so  that  in  general 

we  need  a  four-by-four  matrix  to  describe  the  medium.  However, 
we  usually  assume  a  scalor  behavior  for  the  dielectric,  reducing 
our  four-by-four  matrix  to  two  two-by-two  matrices,  one  tor  TE 
waves,  and  one  for  TM  waves.  Once  we  have  obtained  the  matrix 
coefficients  for  y  and  the  problem  is  essentially  solved. 

3 . 2  Relationships  Between  Two  Different  Approaches 

Since  the  medium  is  periodic,  the  dispersion  curve  is  given 
by  an  expression  of  the  form^^ 

cos  kd  =  Tr  [A]  ,  (3.1) 

where  A  is  a  transfer  matrix.  The  elements  of  A  will  be  func¬ 
tions  of  the  frequency  of  the  incident  wave,  to,  and  the  angle  of 
incidence  of  the  incoming  wave  with  respect  to  the  filter. 
Alternatively  we  may  find  the  input  impedance  and  charac¬ 
teristic  impedance,  Zq  (or  their  corresponding  admittances)  from 
the  elements  b^j  of  A,  from  relations  of  the  form,'^'^ 


and 


"c  =  '^1  -  ‘’22  i  ^  ^  ^  ^22>  .  (3.3) 

In  Equation  (3,2)  Zj^  is  the  load  impedance.  We  see  from  Equation 
(3.2)  that  the  impedance  concept  may  be  used  to  describe  the 
effect  of  terminations,  junctions,  and  local  perturbations. 

However,  these  quantities  can  also  be  described  by  the  difference 
equation  approach.  The  impedance  concept  is  most  useful  when  the 
dielectric  is  a  scalar  (resulting  in  two-by-two  A  matrices)  and 
the  problem  is  one  dimensional.  More  general  problems  can  be 
expressed  in  terms  of  Green's  functions,  which  are  closely 
related  to  A,  and 

From  Equations  (3.1)  (3.3)  we  see  that  we  may  obtain  the  disper 

sion  relation  in  terms  of  the  system  impedance.  If  the  system  is 
viewed  as  a  one-dimensional  transmission  line,  we  have  impedances 


zii,  Z22t  zy2  (see  Figure  3.2),  with^^ 

ail  =  Z22  /  ^12  »  (3.4) 

322  =  ^11  /  ^12  '  (3.5) 

2 

312  (^12  “  2iiZi2)/212  '  (3.6) 

a2i=‘”l/zi2'  (3.7) 

and 

cosh  kd  =  (z]^]^  +  Z22^/2z]^2  •  (3.8) 

Zc  =  Y  (zil  “  "^22^  ^12  sinh  kd  .  (3.9) 


The  reflection  and  transmission  coefficients  may  be  obtained 
directly  from  A,  or  from  We  have‘*^»^0 


Ri 


321  /  322 


(3.10) 


R2  =  ai2  /  322  ' 


(3.11) 


and 


T12  ■  '^21  ~  1/^22  » 

for  reflection  on  the  left,  ,  reflection  on  the  right,  R2 ,  and 
transmission  ”^21.  Of  course,  we  will  also  have^*^ 


R 


z  -z  . 

O  1 

z  +z . 

O  1 


f 


if  the  filter  is  approached  via  an  impedance  Zq. 


(3.12) 


The  dispersion  curve  can  also  be  found  from  a  knowledge  of 
either  8  or  z.  This  is  especially  true  if  A  is  symmetric,  in 
which  case  a]^]^  =  a22»  and/or  an  absorption  may  be  neglected, 
which  implies  that  det[A]  =  1.  Since  the  stop  and  pass  band 
edges  occur  at  cos  kd  =  1 ,  in  the  case  in  which  absorption  may  be 
neglected,  the  dispersion  curve  gives  the  filtering  behavior  for 
a  given  form  for  the  modulation  profile.  This  provides  an 
understanding  of  the  inverted  problem,  the  specifications 
required  on  the  modulation  potential  to  obtain  a  set  of  predeter¬ 
mined  band  edges.  Similar  comments  apply  to  the  reflection  and 
transmission  curves,  provided  again  that  we  have  expressions  in 
terms  of  functions,  rather  than  constants.  The  required 
expressions  will  become  unwieldy  if  the  solution  for  if/  within  the 
unit  cell  of  length  d  cannot  be  expressed  in  a  compact  form. 
Nevertheless,  we  know  enough  from  earlier  work  to  say  a  good  deal 
about  the  behavior  of  the  dispersion  curve  resulting  from  a  given 
modulation  profile. 

For  example,  Yeh,  Yariv,  and  Hong^^  have  studied  light  propa 
gation  in  a  filter  with  a  step  function  ( Kron ig-Penney )  modula¬ 
tion.  The  results  of  their  calculations  are  shown  in  Figures  3.3 
to  3.7.  In  Figures  3.3  and  3.4  the  pass  and  stop  bands  are 
plotted  as  a  function  of  the  angle  of  incidence,  6,  with 


6 


k  sin  e 


(3.13) 

The  case  studied  in  Reference  51  dealt  only  with  small  modu¬ 
lations,  so  the  stop  bands  are  narrow  compared  to  the  pass  bands. 
We  see  that  at  low  angles  of  incidence  the  band  edges  are  insen¬ 
sitive  to  the  angle  of  incidence.  This  insensitivity  extends  to 
larger  angles  as  the  depth  of  modulation  increases.  From  Figures 
3.6  and  3.7  we  see  that  the  reflectivity  curve  corresponds  closely 
to  the  locations  of  the  stop  bands  up  to  fairly  high  angles  of 
incidence  (in  this  case  40").  In  Reference  51  it  is  shown  that 
the  dispersion  curve  (Figure  3.5)  dominates  the  reflectivity 
curve  at  low  to  moderate  angles  of  incidence.  For  TM  waves  the 
stop  bands  narrow  as  they  approach  the  Bragg  angle,  disappearing 
altogether  at  the  Bragg  angle  (Figure  3.7).  We  note  that  Figure 
3.5  shows  how  attenuation  in  the  stop  bands  may  be  incorporated 
into  a  dispersion  curve. 

3 . 3  The  Transfer  Matrix 

Once  we  have  a  solution  for  )(),  we  may  obtain  the  matrix  ele¬ 
ments  of  B  immediately .  52  second  order  PDE  will  have  two  inde¬ 
pendent  solutions.  Let  these  be  and  1^21 


=  Cl 


C2  't>2 


3.14) 


Then 


Cl  !(>i '  +  ^2^12' 


3.15) 


•-  -v*  < 


w (in  units  of 

A 


(a) 


Nfcavcs  at  V 


Refleciivity  spccinim  of  a  1 5-p«nod  Bragg  reflector  for  ( a )  TE  waves  and  t  M  TM 
anous  angle  of  incidence  The  indices  of  refraction  of  the  layers  are  taken  as  n 


}  4 


r  -  thicknesses  are  such  that  a  -  ;  a  The  Brewster  angle  is 


^  .  Reflectivity  as  a  function  of  incidence  angle  for 
wave  in  a  square-well  stack. 


TE 


0.65' 


cos  kd  =  [(}ij^(x)  )  “  ‘{i 


2(x)  (t>'^(x'  ) 


+  i{>^(x)  ())^(x' )  -  4)J^(x)  (|)2(x' )  ]/2D 


(3.27) 


The  expression  of  the  matrix  elements  in  terms  of  x  and  x' 
■  lighlights  their  connection  with  Green’s  functions.  Often  we 
can  reduce  these  expressions  to  expressions  in  terms  of  one  or 
more  distances  between  x  and  x*. 

For  example,  for  a  simple  step  function  modulation  with 
characteristic  distances  a  and  b. 


a  +  b  =  d 


(3.28 


we  have49,53,54 


B(x,x' 


G(a,b) 


H(a,b) 


G'(a,b)  H'(a,b) 


(3.29) 


3 . 4  Constant  Index 

We  now  illustrate  how  this  may  be  done,  using  as  an  example 
the  solution  to  the  equations 


d2ii)/dx2  +  k2||y  =0 


(3.30) 


The  solutions  are 


4>1  =  cos  kx 


<^2  ~  sin  kx 


(3.31) 


(3.32) 


We  let 


x* 


x-a 


(3.33) 


and  find 

D  =  klcos^kx  +  sin^kx]  =  k  ,  (3.34) 

bi2  =  [sin  kx  cos  k(x-a)  -  cos  ka  sin  k(x-a))/k 

=k~^sinka,  (3.35) 

bii  =  [kcos  kx  cos  k(x-a)  +  ksin  kx  sin  k(x-a)l/k 

=  cos  ka  (3.36) 

b2i  =  [-k^sin  kx  cos  k(x-a)  +  k^  cos  kx  sin  k(x-a)]/k 

=  -ksin  ka  (3.37) 

b22  =  [kcos  kx  cos  k(x-a)  +  ksin  kx  sin  k(x-a)l/k 

=  cos  ka  .  (3.38) 

With  the  identifications 

G( a )  =  cos  ka  ,  (3.39) 

H(a)  =  k“l  sin  ka  ,  (3.40) 

we  obtain  the  form  of  Equation  (3.29). 


Expressions  for  the  transfer  matrix  B  in  literature  give  the 
functions  G  and  H  of  Equation  (3.29)  in  terms  of  the  unit  cell 
distance  d: 

G(a,b)  =  G  (a  +  b)  =  G(d)  ,  (3.41) 

H(a,b)  =  H  (a  +  b)  =  H(d)  .  (3.42) 

A  comparison  of  Equations  (3.29)  and  (3.16-3.20)  raises  the 
question  as  to  whether  the  matrix  element  in  Equations 
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I  ^  T I  ^  ■  ■: w  • 


( 3 . 17 ) - ( 3 . 20  )  can  be  expressed  in  terms  of  d.  We  are  also 
iterested  in  the  connection,  if  any,  between  the  G  and  H  func- 
ions  or  equivalently,  between  b^i  and  b22*  such  a  connection 

exists,  we  may  solve  for  the  B  matrix  elements,  given  the  disper¬ 
sion  relation.  This  would  be  a  great  help,  since  we  can  often 
obtain  the  form  for  the  dispersion  relation  directly,  without 
first  finding  the  form  for  B,  as  we  will  show  in  Section  5.  This 
could  be  one  step  in  an  inversion  process. 

Equations  ( 3 . 22 ) -(  3 . 26  )  suggest  a  means  of  answering  both  of  these 
questions.  We  may  rewrite  the  relations.  Equations  ( 3 . 17 ) -( 3 . 20 ) 
in  terms  of  the  difference  between  x  and  x' : 


ii 

il 


►  "  »•  , 


(3.43) 


Usually  we  will  choose  r  =  d.  From  the  definition. 
Equation  (3.43),  we  obtain 


3r/3x  =  1  ,  3r/3x'  =  -1 


(3.44) 


We  write 


■  , 


WL 


and  see  that 


d(|)/dr  =  (}i'  , 


(3.45) 


li  =  ^  il  =  A' 
3x  dr  3x  ^ 


(3.46) 


and 


3 (t>  _  dA  3r 

3x '  ^  TsT" 


(3.47) 


Also, 


2  2 
d^A  3'^r 


3^  A 

3x3x'  ^^2  3x3x' 


(3.48) 


V‘  •«*  ’ 


m 
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It  seems  then,  that  we  may  write 


i.e.,  G  ~  H'  .  (3.50) 

From  Equation  (3.49),  it  appears  that  we  may  write  Equation 
(3.1)  or  Equation  (3.27)  as 

2  b'  =  cos  kd  .  (3.51) 

Consequently,  if  we  have  a  dispersion  relation  of  the  form 

2  g(d)  =  cos  kd  ,  (3.52) 

we  may  write 

b  =  /g  ,  (3.53) 


We  see  then  that  the  central  problem  is  to  solve  the  PDE 
equation  for  ij^.  We  now  discuss  means  of  obtaining  this  solution 
further,  and  give  pertinent  examples  in  which  the  solutions  have 
been  obtained. 
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SOLUTION  OF  THE  WAVE  EQUATION 


We  will  make  some  simplifications  in  our  discussion,  mostly 
for  the  sake  of  clarity.  We  will  confine  most  of  our  remarks  to 
the  one-dimensional  problem.  We  will  also  deal  in  general  with 
real  refraction  indices  and  infinitely  long  systems.  All  of 
these  simpl i f icat ions  can  easily  be  extended  to  treat  a  more 
general  case.  We  will  also  confine  our  discussion  to  systems 
that  are  strictly  periodic  with  period  d.  This  is  an  important 
boundary  condition.  From  the  discussion  that  follows,  we  will 
see  that  it  will  not  seriously  limit  the  performance  of  our 
f ilter . 

We  seek  to  solve  equations  of  the  form^ 


d-^U  ^  ,,2  /  2  ,  .  2  dine  dU 

— 2  -  s  ]U  = 

dx 


dx  dx 


and 


li: 

dx' 


^  ,  ,,2  ,2  _  _2.„  _  din(w-s2/e)  dV 

+  k^  [n  (X)  -  s  JV - ^ ^ 


for  TM  waves  and 


^  +  k^  (n^(x)  -  s'‘)U  = 
dx^  ” 


2 1,,  ^  diny  dU 
dx  dx 


and 


(4.1) 


(4.2) 


(4.3) 


jjv  (n^x)  -  dV  , 

,  2  o  '•  •'  dx  dx 

dx 

for  TE  waves.  Even  if  y  is  a  constant,  e  will  not  be.  If  e 
varies  slowly  we  may  make  a  WBKB  approximat ion^ f 3  g^d  set  the 
right  side  of  these  equations  equal  to  zero.  In  general, 
this  approximation  will  not  be  valid.  It  can  serve  as  a  good 
starting  point  for  perturbation  theory,  it  makes  contact  with  a 
large  body  of  literature,  and  it  is  instructive.  Therefore  we 
will  consider  solutions  to  the  equation 


■  ^  *,  •«*  ■»*  ’•*  *  *  • 
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without  loss  of  generality,  we  may  simplify  our  discussion  by 
setting  s^  =  0  (normal  incidence). 

In  mechanics  we  have'' 

n(x)  ~  /(E  -  V(x) )/E  (4.6) 

where  E  is  the  energy  and  V(x)  is  the  potential.  In  the  par¬ 
ticular  case  of  quantum  mechanics,  we  have 

k  =  p/h  , 

and  an  equation  of  the  form  Equation  (4.5)  which  is 
Schroedinger ' s  equation.  A  great  deal  of  work  has  been  done  on 
periodic  systems  using  this  equation  in  the  band  theory  of 
solids.  n  corresponds  to  the  impedance  z  in  transmission  line 
problems.  In  acoustics,  we  have  n  ~  pv,  for  velocity  v  and  den¬ 
sity  p.5 


4.1  Piecewise  Constant  Modulation 


The  step  function  modulation  is  one  of  the  most  thoroughly 
analyzed  of  any  periodic  response  function. It  is  easy  to 
construct  experimentally,  can  be  solved  exactly,  and  has  cosines 
and  sines  for  eigenfunctions,  so  that  products  of  eigenfunctions 
can  often  be  reduced  to  a  compact  form.  Since  the  system  is 
piecewise  solvable,  we  may  multiply  matrices  of  the  form 


cos  k  n  .  «  . 

oil 


n.  ^  sin  k  n.il. 
1  oil 


-n .  sin  k  n . £ . 
1  oil 


cos  k^n.i. 


(4.7) 


for  each  layer  of  thickness  and  index  n^.  The  most  studied 
case  is  the  double  layer  (one  atom)  ,  with  periods  a  and  b  and 


indices  and  n^.  After  multiplying  matrices,  we  obtain  the 
dispersion  relation  from  Equation  (3.1); 

cos  kd  =  cos  k^n,a  cos  k^n.b 
o  a  o  b 

1  "^a  ^b 

-  4-  f—  +  —1  sin  k^n,a  sin  k^n.b  .  (4.8) 

2  ''n.  n  •'  o  a  o  b 

b  a 

Here  d  =  a  +  b.  We  see  that  when  n^  +  n^  =  n,  we  find 

cos  kd  =  cos  kQn(a  +  b)  ,  (4.9) 

as  we  should.  When  n^  <<  n^  and  a  ~  b  we  find 

cos  kd  =  cos  kQHj^b  -  -j  sin  k^n^b  .  (4.10) 

A  problem  of  more  widespread  interest  occurs  when 
nbb  <<  nga,  but  n^  >>  n^.  Then  we  have 

cos  kd  s  cos  kQn^a  -  (nb/n^)  kbb  sin  konga,  (4.11) 

where  kb  =  ko'^b*  This  system  was  studied  by  Mason  in  regard  to 

acoustic  filters  as  early  as  1927, ^  but  is  best  known  as  the 
Kronig-Penney  problem,  after  these  authors  published  their  inde¬ 
pendent  derivation  of  the  solution  applied  to  a  model  of  energy 
bands  in  solids  in  1931.10  Usually  a  variable  ng  is  defined  in 
such  a  way  that 


n  a  = 
e 


^(nb/na)kbb 


(4.12) 


Using  Equation  (4.6), 
V(x)  =  0  over  region  a,  and 

so  that  i  nbkQb/na  Ub/k 

this  model  is  commonly  used 

If  we  make  the  association 


with  V(x)  =  -U  >>  E  over  region  b, 

2  2 
E  =  kQ  we  have  n^  =  1  and  nb  =  U/E 

Q,  This  is  commonly  written  as  U'a. 

for  energy  band  calculations. 

_  1  ' 

~  -j-  kobnb  =  nb  =  Zb»  we  make 
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contact  with  the  problem  of  a  transmission  line  with  unloaded  line 


impedance  that  is  loaded  at  distances  a  apart  by  discrete 
impedances  z^: 

cos  kd  =  cos  kon^a  -  sin  k^naS  .  (4.13) 

For  a  String  or  rod  loaded  by  point  masses,  Z5  =  mw. ^  If  we  now 
consider  only  long  wavelength  waves,  we  may  write: 

1  -  cos  kd  =  (n^'/n^)  k^nga  ,  (4.14) 

or 

2  sin^  kd)  =  kon^a  ,  (4.15) 

making  contact  with  the  literature  on  lattice  dynamics. 

The  double  periodic  square  well  potential  has  been  studied 
by  Dubrovskii  and  Pogorelsk i i ; the  resulting  expression  for  the 
dispersion  relation  is  complicated. 

4 . 2  Sinusoidal  Modulation 

Another  problem  that  can  be  solved  exactly  is  that  in  which 
the  potential  is  sinusoidal.  This  was  first  treated  by  Morse  in 
1930.29  Morse  solved  the  problem  in  three  dimensions  as  well  as 
in  one  dimension.  We  write  e(x)  as 

e(x)  =  Cq  ~  ejcos  ( 2TTx/a )  .  (4.16) 

and  letting 

a  =  ^  “  *5  *^0^1  '  e  =  (4.17) 

we  may  write  Equation  (4.5)  (with  s  =  0)  as 

d^ij,  /  dx2  =  [a  -  2qcos  ( 2nx/a )  ]  .  (4.18) 

This  is  the  standard  form  of  Mathieu's  equation. 20  The  solutions 
are  circular  elliptic  functions.  There  are  four  basic  types  of 


these  functions; 


ce2n(z»q)  =  1 
r 

^2r  • 

(4.19) 

ce2n+i(z»q)  =  1 
r 

^2r+l  cos(2r+l)z  , 

(4.20) 

se2n+l(2'^^  =  1 

r 

®2r+l  sin(2r+l)z  , 

(4.21 ) 

®®2n+2<=''^>  =  1 

®2r+z  sin(2r+2)z  , 

(4.22) 

where  z  =  irx/a  and  should  not  be  confused  with  the  impedence. 

The  functions  ce^  and  se^,  go  over  to  cosine  and  sine  functions 
when  is  small  compared  to  cq*  Some  of  them  are  shown  in 
Figure  4.1.  The  Mathieu  functions  may  also  be  expanded  in  terms 
of  cylindrical  Bessel  functions  and  Hankel  functions,  and  in  the 
case  of  large  q,  may  be  approximated  by  parabolic  cylinder 
functions. We  will  discuss  various  expansions  of  the  Mathieu 
functions  in  Section  6. 

Levine  has  examined  the  effect  of  various  terminations  on 
the  stop  and  pass  bands  pertaining  to  the  Mathieu  problem. 
Perturbations  were  studied  by  Slater, 32  who  considered  a  dielectric 
of  the  form  Cq  +  eicos  2z  +  e2COS  2pz,  as  well  as  eq  +  e^cos  2z 
+  Slater  considered  two-  and  three-dimensional  per¬ 

turbation  problems  as  well  as  the  one-dimensional  problem. 

4 . 3  General  Observations 

Since  both  the  Kronig-Penney  and  the  Mathieu  problems  are 
exactly  solvable,  comparisons  have  been  made  of  the  stop  and  pass 
bands  created  by  these  modulation  profiles.  The  result  depends 
on  the  relative  size  of  the  distances  a  and  b  in  the 
Kronig-Penney  profile.  When  a  =  b,  the  two  profiles  give  very 


FUNCTION 


ceo(z,q) 


se| (z, q) 


ce|  (z,q) 


se2(z,q) 


ce2(z,q) 


I 


sin  z 


cos  z 


sin  2z 


cos  2z 


se3(z,q) 


sin  3z 


similar  results.  The  width  of  the  pass  bands  relative  to  the 
stop  bands  increases  with  the  ratio  of  b/a.33 

As  a  wave  propagates  through  the  filter,  it  is  partly 
reflected  due  to  modulations  in  the  refractive  index. 

If  the  wavelength  is  much  greater  that  the  periodicity  in 

(nkQd<<l  for  a  mean  index  n),  the  wave  will  not  see  the  index 
modulation,  and  the  wave  will  be  passed.  When  the  wavelength  in 
the  medium  nearly  matches  d,  the  wave  will  interact  strongly  with 
the  modulated  profile;  constructive  interference  will  occur  between 
reflections  that  occur  at  distances  d  apart.  This  will  also 
occur  at  wavelengths  which  are  submultiples  of  this  critical 
wavelength — a  series  of  stop  and  pass  bands  will  occur.  Examples 
corresponding  to  Equation  (4.10)  are  given  in  Figures  4.2  and  4.3, 
where  y  '^b“'^a*  width  of  the  stop  bands  is  roughly  propor¬ 

tional  to  the  amplitude  of  the  index  modulation. 

By  our  argument,  we  see  that  successive  stop  bands  arise  in 
part  from  the  corresponding  Fourier  components  of  the  modulation 
profile.  Thus  the  second  stop  band  is  sensitive  to  the  strength 
of  the  second  Fourier  component  of  the  index  modulation,  etc. 
Consequently  the  widths  and  positions  of  the  lowest  bands  are  not 
very  sensitive  to  changes  in  the  modulation  profile,  only  to  its 
periodicity  (length  of  the  unit  cell).  Changes  in  the  shape  of 
the  profile  have  their  greatest  effect  on  the  higher-order  bands. 
For  example,  if  we  approximate  a  sinusoidal  index  modulation  by  a 
series  of  thin  steps,  we  may  expect  very  little  effect  on  the 
location  of  the  edges  of  the  lowest  bands.  At  high  bands 
corresponding  to  wavelengths  on  the  order  of  the  step  thickness, 
we  would  anticipate  relatively  large  changes.  Of  course,  these 
higher  order  bands  may  well  correspond  to  wavelengths  that  lie 
above  the  bandgap  of  the  dielectric  materials.  If  the  step 
thickness  is  100  nm,  then  we  would  expect  strong  changes  in  the 
bands  corresponding  to  energies  of  say  lOeV.  However,  this 
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d=2a 


Energy  bands  for  various  values  of 
ya,  with  d  —  4a 


Figure  11.  Dispersion  curves  for  d  =  4a. 


region  is  not  of  particular  interest.  In  other  words,  the  modu¬ 
lation  depth  and  overall  periodicity  d  are  the  main  parameters 
affecting  the  band  edges  of  interest.  Fine  details  of  the  shape 
of  the  modulation  are  not  very  important.  Bands  corresponding  to 
the  Mathieu  problem  are  shown  in  Figures  4.4  and  4.5.  A  com¬ 
parison  between  the  bands  for  the  Kronig  Penney  and  Mathieu 
problems  is  shown  in  Figures  4.6  and  4.7. 

If  the  system  is  multiple,  periodic  higher-order  bands  at 
corresponding  multiples  will  be  affected.  We  will  discuss  this 
further  in  the  next  section,  but  we  can  see  from  the  argument 
given  above  that  since  a  multiple  periodicity  implies  that  second 
and  higher  order  Fourier  components  of  the  modulation  profile 
will  be  large,  the  stop  bands  (band  gaps)  above  the  fundamental 
will  be  widened. 
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5.  DIFFERENCE  EQUATION  CALCULATIONS 

Often  we  can  convert  a  PDE  into  a  difference  equation  by  an 
appropriate  series  transformation  or  decomposition,  followed  by 
an  approximation  which  allows  us  to  limit  the  coupling  that 
occurs  between  the  basis  states  of  the  decomposition.  We  now 
give  two  examples  of  how  this  can  be  done.  The  first  is  a 
Fourier  transformation;  the  second  is  known  as  the  tight  binding, 
or  LCAO  method.  It  is  the  approximation  which  limits  the  number 
of  connected  basis  states  appearing  in  the  difference  equation 
that  makes  the  difference  equation  method  practical.  consequently, 
this  approach  is  only  useful  when  such  an  approximation  can  be 
made . 


5.1  Fourier  Transformation^*^ 

We  wish  to  solve  the  PDE  of  the  form  in  Equation  (4.5).  We 
make  the  decompositions 


ik  T  X 

=  1  e  ^  » 

k  j  j 

where 


kj  =  k  +  2Trj/d  , 


and 


2 

n 


(X) 


I  g 

I 


iK^x 


r 


(5.1) 


(5.2) 


(5.3) 


where 


K£  =  2Ti£/d  . 


Substituting  Equations  (5.2)  and  (5.3)  in  (4.5),  we  find 
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I[(k  .  .  k^[s^  -  I  g, 


)]  A.(k)e 


(5.4) 


In  principle,  Equation  (5.4)  is  exact,  but  involves  an  infi¬ 
nite  sum  over  the  Fourier  components  g^  for  the  square  of  the 

index,  n2^  need  to  drastically  reduce  the  number  of  these 

2 

components.  Fortunately,  if  n  (x)  is  a  reasonably  smooth 
function,  we  may  terminate  the  series  in  I  after  a  few  terms.  As 
an  example,  we  consider  the  Mathieu  problem,  in  which  we  have 


(5.5) 


Substitution  of  Equation  (5.5)  into  Equation  (5.4)  now  gives 


f  f  k  +  k.]  ^  +  k^s^l  A.  =  .k^[A.  T  ] 

^  ^  o  J  0  J-1  D+1 


(5.6) 


This  is  a  solvable  equation.  Rather  than  pursue  it  here,  we  go 
on  to  derive  a  similar  equation  by  another  means. 

5.2  Tight-Binding  Approximation3 

A  great  deal  of  work  has  been  done  using  difference  equation 
calculations  in  the  band  theory  of  solids  based  on  the  tight- 
binding  approximation.  3-12  xhis  corresponds  to  the  case  in  which 
regions  of  large  refractive  index  are  sufficiently  separated  by 
regions  of  low  refractive  index  that  they  form  optical  wells 
which  are  only  weakly  coupled.  Parallel  calculations  have  been 
made  in  three  dimensions  as  well  as  one  dimension  and  include 
investigation  of  the  effect  of  various  terminations  (surface 
states).  Although  they  are  most  accurate  in  cases  of  very  deep 
modulation  of  the  index  profile,  they  provide  an  instructive 
insight  into  band  behavior,  and  serve  as  good  points  for  com¬ 
parison  with  more  general  calulations. 


Let  us  assume  that  we  have  a  complete  set  of  known  orthonor¬ 
mal  solutions  (|)  „  for  the  equation 


”22  2 

(t  +  k[n{x)-s]ii)  =  0  , 

o*-  o  O'* 

and  that  we  wish  to  solve  Equation  (4.5): 

"  +  k^[n^(x)  -  =  0  , 


(5.7) 


(5.8) 


for  some 


n^(x)  =  n^(x)  +n^(x) 


(5.9) 


We  may  expand  iJj  in  terms  of  some  as' yet  unknown  coefficients 


as  follows: 


'P  =  I  a. 


(5.10) 


Upon  substitution  of  Equation  (5.10)  into  Equation  (5.8)  we  obtain 


2  2  2  2  2 

Takfs  -  sli(i  =yakn4> 

^  n  o'-  o-'  n  ^  n  o  n 

n  n 


(5.11) 


If  we  multiply  this  equation  of  <|in,  and  integrate,  we  will  obtain 
the  results 


[  a  k^[  s^  ■  =  I  3_k^dx  ip  (x)<p  (x)ri^(x) 

‘■no'-  o-'  ‘-no^  ra  n 

n  n 


(5.12) 


Since  the  functions  iti^(x)  are  known,  the  integrals  are  obtainable 
in  principle,  and  we  may  solve  for  the  coefficients  a^^  in  terms 
of  them.  The  result  will  be  a  series  of  coupled  equations.  This 
can  get  very  messy  if  a  given  initial  solution  is  coupled  by 
n(x)  to  very  many  other  solutions  (fp,.  The  conversion  to  a  dif¬ 
ference  equation  is  only  practical  if  the  coupling  is  only  to  a 
few  other  solutions.  An  example  in  which  this  can  occur  is  when 
the  <p  ^  are  localized  in  a  region  smaller  than  the  (unit  cell) 
distance  d.  We  consider  the  case  in  which  we  have 

/  dx  ^n,(x)<^^(x)n2(x)  =  a  ,  m  =  n 


=  S  .  m  =  n  +  1 


=  0  ,  (m  -  n)  >  1 


(5.13) 


Then  we  will  find 


e(a  +  a  .  )  =  (s  -  oi)a 
n+1  n-1  n 


(5.14) 


with  s2  -  Sq  =  s.  We  see  that  we  have  greatly  simplified  the 
solution  by  our  approximation.  It  has  been  found  useful  to 


assume  that  a^  is  of  the  form 


an  =  y 


(5.15) 


Substitution  of  Equation  (5.15)  into  Equation  (5.14)  yields 


[6X2  -  (s  -  a)X  +  6]  =  0 


(5.16) 


We  now  rewrite  Equation  (5.16)  as 


-  [^]X  +1  =  0. 


(5.17) 


We  note  that  Equation  (5.17)  is  of  the  same  form  as  the  eigenvalue 
equation  for  the  transfer  matrix,  B.  That  is,  if  we  wish  to 
solve  the  equation 


4)(x+d)  =  Xi);(x)  , 


(5.18) 


we  need  to  find  the  roots  of  the  transfer  matrix  as  given  for 
example  in  Equation  (3.28): 


G(d)  -  X 


G'  (d) 


H(d) 

H' (d)  -  X 


(5.19) 


This  gives 


a2  -  [G(d)  +  H’ (d)  ]X  +  1  =  0  , 


(5.20) 


since  det ( 8 1 


1. 


We  know  that  the  dispersion  relation  is 

Tr[B]  =  G  +  H'  =  2coskd  ,  (5.21) 

and  take  advantage  of  this  insight  to  make  the  assumption 


(s  -  a)/B  =  2coskd  .  (5.22) 

Now  if  we  substitute  Equation  (5.22)  into  Equation  (5.17),  we 
will  find 

X  =  coskd  isinkd  =  e- ^  (5.23) 

In  turn,  substitution  of  Equation  (5.23)  into  Equation  (5.15) 
will  bring 

a^^  =  peiknd  +  Q^-iknd  ^  (5.24) 

Since  the  basis  functions  <)>  ^^  were  assumed  to  be  known,  the 
problem  has  now  been  solved  given  the  appropriate  boundary  con¬ 
ditions.  We  see  from  Equation  (5.24)  that  we  have  shown  that  the 
assumption  Equation  (5.22)  implies  that  the  expansion  Equation 
(5.10)  is  in  terms  of  Block  waves.  Usually  the  equivalence  is 
derived  in  the  other  direction:  the  Block  expansion  allowed  due 
to  the  periodicity  of  the  system  implies  the  condition  Equation 
(  5 . 2 1 )  .  1 3  We  note  that  Equations  (5.6)  and  (5.14)  are  of 

similar  form.  This  is  in  part  due  to  our  use  of  a  Block 
expansion. 


5.3  Surface  States 


We 

through 

periods 


can  investigate  effects  such  as  terminating  the  filter 
the  boundary  conditions.  We  let  the  filter  consist  of  N 
of  length  d,  so  that  the  filter  has  a  length 


59 


L 


Nd 


(5.25) 


We  then  impose  the  boundary  condition  that  be  zero  at  L: 

=  0  .  (5.26) 

Together  with  Equation  (5.24),  Equation  (5.26)  implies  that 

a^  =  Asin(N  -  n)l<d  /  sinNkd  .  (5.27) 

Since  A  /  sinNkd  is  simply  a  constant  normalizing  factor,  we 
drop  it  in  favor  of 


We  obtain  our  other 
making  the  substitution 

eai  =  (s  -  a ' )aQ  (5.29) 

for  Equation  (5.14).  Then  Equations  (5.28)  and  (5.29)  give 


=  sin ( N  -  n) kd  .  (5.28) 

boundary  condition  at  the  termination  by 


(s  -  c')sinNkd  =  6sin(n  -  l)kd  .  (5.30) 

From  Equations  (5.30)  and  (5.22)  we  get 

(a  -  a')/B  +  2coskd  =  coskd  -  cotNkd  sinkd  ,  (5.31) 

or,  with 

^  =  (a  -  a ' )  /  6  ,  (5.32) 

C  =  -  sin(N  +  l)kd  /  sinNkd  .  (5.33) 


The  intersection  of  this  family  of  curves  with  the  curves 
;(kd)  =  constant  yield  solutions  for  c  >  1.  One  solution, 
corresponding  to  a^,  the  "surface  state,"  will  give  a  complex 
value  of  k.  This  wave  will  not  propagate  through  the  filter. 
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The  difference  equation  tack  is  a  useful  means  of  studying 
the  effects  of  multiple  periodicity.  The  conversion  to  a 
difference  equation  has  been  used  extensively  for  this  purpose  in 
band  theory  and  lattice  dynamics.  The  published  literature  is  a 
valuable  guide  for  parallel  calculations  which  we  might  wish  to 
make.  However,  the  results  that  have  been  obtained  in  these 
fields  are  not  readily  used  in  optics,  since  energy  band  results 
seek  to  solve  for  the  eigenvalues  (or  energy)  of  the  Schroedinger 
equation,  which  is  equivalent  to  solving  for  s  in  Equation  (5.12) 

or  s  in  Equation  (5.14).  We  can  easily  take  s  =  0,  and  wish  to  find 
the  relation  between  k^  and  k.  We  are  also  interested  in  an  algorithm 
that  is  generally  useful  beyond  the  nearly  free  and  tight  binding 
approximations.  It  seems  that  this  can  be  done.  Let 


<l>  '  (x) 


^  (x) 


(5.34) 


Then 


<ti  (x+d)  =  B(d)i|)  (x) 


lii(x-d)  = 


Hence 


Wi  th 


B~^(d);i;(x) 


i|;(x-fd)  +  ij;(x-d)  =  [B(d)  +  B  ^(d)]t(((x) 


(5.35) 


(5.36) 


(5.37) 


k'-j'  •y'J 


ij/  (  x+d  )  =  e  !(/  ( X ) 


(5.38) 


/  ~kd  ,  , 

4)  (  x-d  )  =  e  t  (  X )  > 


(5.39) 


we  find 


B(d)  +  B  (d)  =  2cosh(kd) 


(5.40) 


From  Equation  (3.29)  (and  the  discussion  following  it  in  Section  3) 


we  have 


G'  H' 


(5.41) 


The  inverse  of  B  is 


’  H' 

5-1  = 

_-G' 


(5.42) 


Therefore  we  may  write  Equation  (5.40)  as 


G+H' 


=  2  cosh(kd)  , 


(5.43) 


(G+H'  ) 


0  1 


=  2cosh(kd) 


(5.44) 


It  seems  that  the  difference  method  may  be  applied  to  a  general 
problem,  if  the  boundary  conditions  are  suitably  chosen  with 
respect  to  the  unit  cell.  This  would  be  very  useful  for  the 
analysis  of  a  multiply  periodic  medium.  Consider  the  system 


4*k  =  ex 


tl*-  =  Ot  X 


(5.45) 


(5.46) 


where  the  subscripts  refer  to  distances  with  respect  to  the  argu¬ 
ment  of  the  r.h.s.  We  find 


/.  — I'.'*' 

'i'a  =  “  Jx 


Xa  =  a  4- 


Xg  =  0  I* 


Therefore 


X;,  +  Xh  =  (a  +  0  )’!' 


A.a  ■  .  K  ;r  • 

From  Equations  (5.47)  and  (5.50)  we  get 

(S  +  a  (a  +  0  =  (e  +  a  (Xg  Xg) 

=  '('h+fl  +  ^  +  I 


We  write  Equation  (5.51)  as 


With 


[(6  +  a  ^)(o  +6  -  2]i);  =  +  'I'g 


4'^  =  X4; 


4'g  =  X  ^4> 


we  have 


(5.47) 


(5.48) 


(5.49) 


(5.50) 


(5.51) 


(5.52) 


(5.53) 


(5.54) 


*•* 

liir'l 
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(B  +  a  )(a  +6  )  -2=A  +A 


(5.55) 


or 


Sot 


X  + 


We  let 


M  =  6a 


and 


Then  we  have 


Tr[M]  =  coshkd 


(5.56) 

(5.57) 

(5.58) 

(5.59) 
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6. 


PERTURBATION  EXPANSIONS 


We  now  look  at  some  of  the  ways  in  which  we  can  find 
approximate  solutions  for  iji  and  for  the  dispersion  relation  when 
the  index  modulation  has  a  small  departure  about  a  profile  for 
which  we  can  solve. 

6 . 1  Pedestrian  Expansion 

If  we  know  that  our  modulation  profile  is  similar  to  one  for 
which  we  have  a  closed  form  solution  we  may  make  expansions  about 
both  the  known  solution  and  its  corresponding  modulation  profile. 
We  let  the  known  solution  be  i(;  ^ ,  and  the  unknown  solution  be  ip . 

We  then  write 


^  ^  ^  ^1*3  '  1  •  (6.1) 

As  5  increases,  we  will  need  to  include  more  terms  in  the  expan¬ 
sion  for  tp ,  Similarly,  if  the  dielectric  constant  for  the  known 

2 

case  goes  like  n2(xl),  we  write 

n^(x)  =  nj(x)  +  ^  (6.2) 

for  n2(x)  in  the  unknown  case.  Then,  upon  substituting  Equations 
(6.1)  and  (6.2)  into  Equation  (4.5)  with  s^  =  Q; 

d^i|</dx^  +  k^n^{x)i(<  =  0  ,  (6.3) 

we  obtain 


k^(n^  +  6  n^  +  6^n^)(il;,  +  ~  ^  •  (6.4) 


We  can  separate  terms  of  a  given  order  and  arrive  at  a  series  of 
equations : 

ij,  J  1  =  0'  (6-5) 

’^2  ~  '^2’*' 1  '  (6.6) 

and 

Equation  (6.5)  is  our  solvable  homogeneous  equation  with  a 
known  solution.  Equations  (6.6)  and  (6.7)  are  similar  to 
Equation  (6.5)  but  are  inhomogeneous.  We  can  think  of  Equation 
(6.6)  as  a  problem  in  which  the  solution  to  Equation  (6.5)  drives 
an  oscillator  of  the  same  form  as  that  of  Equation  (6.5). 

We  can  learn  something  about  the  dispersion  relation  (k^ 
versus  k)  if  the  amplitude  of  the  index  modulation  is  small  by 
making  another  expansion,  similar  to  that  of  Equation  (5.1).  We 
write 


ik  .  X 

'J'l,  (x)  =  I  A.  (k)e  ^ 


J 


J 


(6.8) 


with 


k  j  =  k  +  2ii  j/d 


(6.9) 


6 . 2  Band  Gap  for  Small  Perturbations 


We  assume  a  small  perturbation  in  the  index  modulation  about 
a  value  for  which  we  can  solve.  Then  we  have2f3 


Ij)  (X  =  A  e  +  6  )  A.e 

k  m  Ir  j 

We  begin  by  letting  m  =  0.  We  write  our  equation  for  i|;  as 


(6.10) 


(6.11) 


with 


_  iK  X 

2  r  Z 

n,  =  ig.e  ■ 

^  I  ^ 


(6.12) 


=  2vZ  /d 


(6.13) 


Then  we  substitute  Equations  (6.9),  (6.10),  and  (6.12)  into 
Equation  (6.11)  and  find 


A  (  k^n^  -  k^]e^'^’^  +  <5  -  k'^)A.e  ^ 

o''  O  O  ‘  1“  '•  o  o  v  j 

D  J  J 


V.2_2  .2, 


ik  .  X 


T  X 

+  6 .k^A  y g„ e 
2  o  o^  ^Z 
Z 


+  5,6  _k^y  y  g„  A  .  e 
1  2  oA  j 


i(k^+K  Jx 


=  0. 


(6.14) 


To  first  order  we  let 


Then  to  first  order  we  have  approximately^ 

A  (k^n^  -  k^je'"’' 

0^0  0  ^ 


*  - "?)  ''j  *  ^  '  0 

1 


ik  .  x 


(6.15) 


We  now  operate  on  Equation  (6.15)  with 


rJ*  -m-r  WJ v; y; »v ’>'.' *V »'.' ^T,' n.' jnii>«jr»w 'i^i^ v.- « ■» »/  «.i» -j-^-^TT^-TTy 


ik  X 

,  .  „ikx  ^  -  m 

ip  =  A  e  +  A  e 

o  m 


ik  .  X 


+  6  J"A,e  J 


(6.22) 


Then  proceeding  as  before,  we  find  to  first  order 


A  f  n^k^  -  k'^l  +  <Sk'^A  g  =0 
o'*  o  o  ‘  o  m-'n 


(6.23) 


and 


g  +  ( n^k^  -  ]  A  =0  . 

o  o  -m  '■  o  o  •  -m''  m 


(6.24) 


Hence 


[n^k^  -  k^]  ~  1  g 

^  o  o  ■'  ''  o  o  -m''  o  * 


(6.25 


For  k  ^  k_j^  we  find  approximately 


1=0  ki  /  (nj  .  6g^)  . 

m 


(6.26) 


This  creates  a  band  gap  with  width 


g  g 

Ak  =  — ^  k  =  J  k  . 
..  2  m  2  £  m 

2n  o 

o 


(6.27) 


This  means  that  a  stop  band  is  formed  whose  width  for  the  mth 

2 

band  is  proportional  to  the  mth  component  of  n^^. 
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7.  CANDIDATE  MATERIALS 

The  impetus  for  this  report  was  a  literature  survey  of 
suitable  candidate  materials  for  graded  index  optical  components 
and  devices.  Such  an  undertaking  would  be  a  prodigious  task, 
especially  in  view  of  the  limited  data  available  on  the  optical 
properties  of  materials,  and  the  vast  number  of  possibilities. 
This  is  all  the  more  true  because  of  the  decision  to  direct  the 
bulk  of  the  effort  of  this  work  to  the  analysis  discussed  in  the 
first  six  sections  of  this  report.  Nevertheless,  the  author 
wishes  to  present  his  initial  impressions  of  the  possibilities 
for  suitable  materials,  albeit  from  the  prospective  of  a  novice 
and  an  outsider,  and  to  report  on  a  limited  literature  survey. 
First,  some  general  considerations  and  comments  are  in  order. 

7.1  General  Comments 


The  throughput  of  an  optical  device  and  its  immunity  to 
damage  under  high-power  conditions  will  be  enhanced  if  the 
absorption  is  low  for  a  given  optical  path  length,  i.e.,  if 
ni  <<  n^.  This  mitigates  against  the  use  of  amorphous  materials, 
which  in  practice  tend  to  suffer  from  large  absorption.  This  is 
due  in  part  to  a  very  high  density  of  scattering  centers  arising 
from  local,  microscopic  inhomogeneities,  such  as  the  frequent 
occurrence  of  voids,  and  a  general  tendency  to  a  high  density  of 
impurities.  This  is  especially  likely  if  the  stoichiometry  must 
be  varied  over  a  wide  range.  The  latter  problem  is  no  doubt 
enhanced  by  the  former  difficulties.  Furthermore,  amorphous 
materials  tend  to  possess  wide  regions  over  which  infrared 
radiation  is  absorbed  by  photons.  Since  we  wish  to  vary  the 
refractive  index  continuously,  a  second  implication  of  the  need 
to  minimize  absorption  and  scattering  is  the  high  desirability  of 
using  materials  with  good  mutual  solubility;  we  seek  compounds 
possessing  a  similar  structure  with  as  different  a  refractive 
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index  as  possible.  The  third  corollary  is  the  requirement  that 
the  chosen  compounds  be  amenable  to  manufacture  with  a  low  con¬ 
centration  of  defects. 

The  chosen  materials  should  be  rugged,  that  is,  they  should 
be  chemically  resistant,  especially  to  moisture,  should  have  good 
mechanical  strength,  and  a  high  phase  transition  temperature. 
Their  refractive  index  should  not  vary  rapidly  with  temperature. 
The  desirability  of  minimizing  short-range  disorder  suggests  that 
the  solid  solute  materials  be  chosen  with  nearly  equal  bond 
lengths  as  well  as  identical  structures.  This  will  be  especially 
true  if  the  creation  of  a  single  crystal  device  is  attempted.  A 
basic  premise  of  the  following  presentation  is  that  the  bonding 
lengths  will  usually  be  dominated  by  the  nonmetals  in  the  com¬ 
pounds,  rather  than  the  metals. 

This  can  be  seen  from  Table  1,  which  gives  the  radii  of 
several  atoms  of  interest.!  We  note  that  covalent  bond  lengths 
are  especially  likely  to  be  dominated  by  the  identity  of  the  non- 
metal  atoms.  In  these  cases  the  small  metal  atoms  often  fit  into 
voids  created  by  nearly  close-packed  nonmetal  atoms.  Thus  we  may 
expect  to  vary  physical  properties  without  significantly 
affecting  bond  lengths  by  substituting  different  metal  atoms  with 
a  given  nonmetal  atom.  A  good  example  of  this  occurs  for  com¬ 
pounds  of  transition  metal  atoms  with  chalcogenide  atoms. 

Assuming  that  this  hypothesis  is  true,  we  should  seek  compounds 
with  identical  structures,  different  cations,  and  anions  which 
are  either  identical  or  of  nearly  the  same  size. 

To  ensure  structural  compatibility  in  a  solution,  we  will 
also  restrict  ourselves  to  groups  of  atoms  which  form  compounds 
with  the  same  structure  as  well  as  closely  matched  lattice 
parameters  when  all  possible  combinations  of  atoms  are  con¬ 
sidered.  Therefore,  we  will  organize  our  discussion  in  terms  of 
crystal  structure  in  general,  and  in  terms  of  the  chemical 


properties  of  the  anions  in  particular.  Since  we  did  not 
consider  dispersion  in  the  earlier  sections  we  will  not  give  it 
much  attention  below.  Also,  we  will  concentrate  on  linear 
refractive  indices.  Finally,  we  will  concentrate  on  simple  com¬ 
pounds,  as  they  should  have  fewer  absorption  bands  in  the 
infrared . 

The  optical  behavior  of  a  compound  is  strongly  affected  by 
the  chemical  nature  of  its  bonds,  in  an  analogous  manner  to  many 
of  its  other  properties.  It  appears  to  the  author  that  the  ionic 
materials  often  have  the  desirable  attributes  of  low  optical 
absorption,  low  defect  concentrations,  and  high  melting  points, 
as  well  as  wide  band  gaps.  In  many  instances,  especially  in  the 
case  cf  the  halogens,  they  are  subject  to  chemical  attack,  par¬ 
ticularly  by  water,  and  often  they  are  brittle.  However,  a  more 
se.  ious  difficulty  is  that  ionic  compounds  tend  to  have  low 
indices  of  refraction,  and  a  range  of  refractive  indices 
extending  to  high  values  is  often  desirable. 

7 . 2  Halides 

The  fluoride  compounds  form  a  group  of  highly  ionic 
materials  that  have  well  developed  and  well  understood  optical 
properties.  They  have  very  low  absorption  coefficients. 

However,  their  refractive  indices  are  uniformly  low  (nj.  ranges 
from  1.25  to  1.53)  seriously  restricting  their  versatility.  Some 
of  the  pertinent  properties  of  several  fluoride  compounds  are 
given  in  Table  2.  In  the  tables,  a  is  the  lattice  constant  (all 
lattice  constants  in  the  tables  are  given  in  angstroms),  is 
the  high-frequency  dielectric  constant,  is  the  static 
dielectric  constant,  and  w  and  are  respectively  the  trans¬ 
verse  and  longitudinal  optical  photon  frequencies  at  the  gamma 
point  in  cm“l.  Gisin^  found  that  films  of  SrF2  showed  improved 
performance  with  increasing  substrate  temperature  during 
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TABLE  2 


FLUORITE  STRUCTURE 


Compound 

a 

n(ym) 

.5 

2 

CdF2 

5.338 

1. 

CaF2 

5.463 

1.437 

1.43-1.46 

HgF2 

5.54 

EuF2 

5.796 

SrF2 

5.800 

1.25-1.36 

PbF2 

5.9  27 

BaF2 

6.200 

1.478 

1. 38-1.53 

SrCl2 

6.977 

BaCl2 

7.34 

5 

1.53 

1.424 


1.451 


yl,-f 


1.458 


fomation.  This  was  attributed  to  a  decrease  in  porosity. 

Gisin  states  that  the  optical  constants  of  films  of  SrF2  depend 
on  the  layer  thickness  due  to  nonuniform  stoichiometry.  Hall3 
found  that  the  index,  absortion,  and  crystallinity  of  MgF2 
increased  with  age.  Hall  attributes  the  increased  absorption  to 
the  creation  of  vacancies  with  the  onset  of  crystallization. 
Barnes  and  Gettemy^  reported  on  the  optical  properties  of  LiYF4/ 
which  has  a  very  low  temperature  coefficient. 

The  CsCl  materials  also  have  low  dielectric  constants.  Data 
on  these  compounds  are  given  in  Table  3,  much  of  which  is  taken 
from  Brown ,  5  Burstein  et  al.,^  Lucovsky,  et  al.,"^  and  to  a 
lesser  extent,  Rao.^  (These  authors  are  also  the  main  basis  for 
Tables  4,  15,  and  24.)  The  best  bet  for  this  group  of  materials 
appears  to  be  either  TICI/NH4CI  or  TlCl/TlBr.  There  is  a  wider 
selection  of  NaCl  structure  materials  with  a  fairly  wide  range  of 
dielectric  constants  (see  Table  4).  Unfortunately,  the  alkalai 
halides  can  develop  color  centers. 5  This  may  present  a  problem 
for  work  in  the  visible,  but  the  Agj^Na^.^Br  and  Ag^Na^_j^Cl 
systems  appear  to  be  attractive  for  work  in  the  infrared.  In 
most  ionic  compounds,  we  expect  only  a  10%  change  in  the 
dielectric  constant  for  two  materials  with  closely  matched  lat¬ 
tice  parameters,  but  in  these  cases  we  find  a  change  of  almost  a 
factor  of  2. 


Compounds  of  the  halides  with  many  of  the  transition  metals 
have  a  layered  structure.  In  these  materials  the  small  cations 
are  trigonally  coordinated  with  the  larger  anions  (see  Figure 
7.1a).  There  are  also  layered  compounds  of  chalcogenide  atoms. 
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TABLE  3 


CsCl  STRUCTURE  HALIDES 


Compound 

a 

^  00 

^  s 

(1)  'p 

“L 

TlCl 

3.84 

4.76 

31 

62 

NH4CL 

3.87 

TLBr 

3.97 

5.34 

CsCl 

4.11 

2.62 

7.20 

99 

165 

Til 

4.20 

29.6 

CsBr 

4  .  29 

2.83 

6.51 

73 

112 

Csl 

4 . 56 

3.02 

5.65 

62 

85 

TABLE  4 


ROCKSALT  STRUCTURE  HALIDES 


Compound 

a 

e® 

Li  F 

4.02 

1.93 

9.27 

Li  H 

4.08 

3.6 

12.9 

Na  F 

4.67 

1.75 

5.3 

Li  Cl 

5.14 

2.75 

11.05 

K  F 

5.33 

1.85 

6.05 

LiBr 

5.49 

3.17 

•  13.2 

AgCl 

5.  55 

4.04 

12.3 

NaCl 

5.64 

2.3 

5.62 

Rb  F 

5.63 

1.93 

6.5 

AgBr 

5.77 

4.62 

13.1 

NaBr 

5.94 

2.61 

5.99 

Lil 

6.00 

11.03 

CsF 

6,01 

KCL 

6.28 

2.15 

4.68 

Nal 

6,46 

2. 96 

6.60 

RbCl 

6.54 

2.18 

4.9 

K  Br 

6.59 

2.35 

4.78 

Rb  Br 

6.87 

2.34 

4.9 

r»T»i 


the  sandwich.  Only  a  weak  Van  dec  Walls  interaction  occurs 
between  the  anions  in  adjacent  sandwiches.  For  this  reason  the 
layered  compounds  are  attractive  in  designing  materials  with 
modulated  physic  properties.  Even  if  adjacent  sandwiches  were 
comprised  of  different  compounds,  rather  than  true  solutions, 
there  should  be  a  low  level  of  interface  states. 

The  arrangement  of  atoms  within  any  given  sheet  is  hexagon¬ 
al.  The  stacking  of  adjacent  sheets  in  layers  that  assume  a 
trigonal  coordination  is  ABC,  so  that  a  plane  taken  through  the 
material  as  shown  in  Figure  7.2  will  contain  an  atom  from  each 
sheet  in  the  structure,  resulting  in  the  representations  given  in 
Figure  7.3  and  7.4.‘^~11  These  figures  depict  the  two  stacking  poly 
types  assumed  by  halide  compounds,  the  Cdl2  and  CdBr2  structures. 
Lattice  constants  for  halide  compounds  with  these  structures  are 
given  in  Tables  5  and  6.  On  the  basis  of  their  lattice 
constants,  Cdl2  type  systems  of  interest  appear  to  be  Ti  ^V]^_j^Cl  2 ; 
bromides  of  Fe-V-Mg-Mn  solutions;  COj^Fej^_2^  1 2 ;  and  iodides  of 
Ti-Ge-Mg-Mn,  Ca-Yb-Tm-Pb,  and  Cd-Zn-Bi-Te  solutions.  Similarly, 
CdBr2  structure  materials  worth  looking  into  are  Znj^Cd  j^_^Br 2 ,  and 
chlorides  of  Ni-Co-Fe-Mg  and  Mu-Zn-Cd  solutions. 

Almost  all  of  the  optical  measurements  on  layered  compounds 
have  been  made  at  short  wavelengths,  with  a  view  toward  eliciting 
information  on  their  band  structures.  Sinha  and  Mukherjee^^ 

Doni  and  Grossol^  have  made  such  measurements  on  Pbl2^  which  has 
many  stacking  polytypes,  and  an  index  of  refraction  of  2.9. 

Yashiro,  et  al.,  reported  on  the  emission  spectra  of  Pbl2f^‘^  and 
Ghita,  et  al  ,  made  el  1  ipsome t r ic  studies  on  it.^-^  References  16 
and  17  report  on  the  phonon  spectrum  of  Pbl2-  Kondo  and 
Matsumoto  made  reflection  measurements  above  4eV  on  CdCl2  and 
CdBr2  crystals,  while  Srivasta  and  Bist  measured  the  infrared 
absorption  of  Cdl2'^^  Anderson  and  Lo^^  give  the  phonon  spectra 
of  several  layered  halides.  Many  trichlorides  also  form  layered 


compounds.  The  infrared  properties  of  TiCi3  were  studied  by 
Erne  is  ,  et .  al .  21 

7 . 3  Oxides 

The  oxides  form  a  natural  group  of  optical  materials,  due  to 
their  wide  band  gaps,  high  melting  points,  and  resistance  to 
chemical  attack.  There  are  far  too  many  of  them  to  be  given  a 
proper  treatment  in  this  paragraph;  consequently  our  discussion 
will  only  be  cursory.  As  elsewhere  in  this  section,  we  will  con¬ 
fine  our  remarks  to  simple  compounds.  The  simplest  system  is  the 
Rocksalt  structure.  There  are  a  large  number  of  oxides  in  this 
system,  with  an  almost  continuous  distribution  of  lattice 
constants,  as  can  be  seen  from  Table  7.  The  most  closely  matched 
systems  are  NixTii_xO,  NbxMgi_xO,  and  PuxPai_xO.  The  latter 
systems  have  a  lattice  match  to  within  0.05%.  Many  of  these  oxides, 
especially  those  of  vanadium  and  titanium,  form  a  series  of 
defect-structure  nonstoich iometr ic  compounds,  Nonstoichiometric 
oxides  are  discussed  at  length  in  Sorensen, 22  and  will  not  be 
dealt  with  here.  We  note  from  the  entrees  on  the  dielectric 
constants  in  Table  7  that  the  oxides  tend  to  have  somewhat  larger 
dielectric  constants  than  the  halides. 

The  Rutile  family  is  a  second  collection  of  familiar  oxides. 
These  materials  are  listed  in  Table  8,  from  which  we  see  that 
the  Ge-Mn-Cr,  Ta-Sn-Nb,  and  WxMoi-x  systems  are  the  best  possibil¬ 
ities  for  Rutile-type  oxides.  As  in  the  case  of  the  dichalcoge- 
nides,  which  we  will  discuss  in  paragraph  7.4,  the  small 
molybdenum  and  tungsten  cations  form  dioxides  with  nearly  equal 
lattice  constants.  Lattice  constants  for  Table  8  and  most  of  the 
following  tables  were  obtained  primarily  from  Wyckof f . 23 , 24  ^ 

similar  but  slightly  more  complex  family  of  oxides  is  given  in 
Table  9,  which  suggests  that  solutions  of  Cr-Ga-Rh-Fe  with  Sb04 , 
of  Nb-Ta  with  Cr04 ,  of  Sb-V  with  Rh04 ,  TaxNbi_xFe04 , 

FexRbi_xNb04 ,  and  FexRhi_xTa04  are  suitable  candidate  materials 
for  continuously  modulated  devices. 
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TABLE  7 


V  and 


ROCKSALT  STRUCTURE  OXIDES 


Compound 

a 

VO 

4.062* 

NiO 

4.  1684t 

TiO 

4.1766* 

NbO 

4.2097 

MgO 

4.2112 

2.95 

CoO 

4.2667t 

TaO 

4.43 

MnO 

4.4448 

4.8 

ZrO 

4.62 

CdO 

4.6953 

5.40 

CaO 

4.8105 

3.33 

YbO 

4.86 

PuO 

4.959 

PaO 

4.961 

SmO 

4.9883 

NpO 

5.01 

EuO 

5.1439 

SrO 

5.1602 

3.46 

Ti 


tend  to  form  a  series  of 


non-stoch iometr ic  oxides 
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TABLE  8 


C^SSITERITE  OR  RUTILE  STRUCTURE  OXIDES 


Compound 

a 

c 

GeO  2 

4  .  395 

2.859 

Mn02 

4  .  396 

2.87  1 

Cr02 

4.41 

2.91 

Ir02 

4.49 

3.  14 

OsO  2 

4.51 

3.11 

RUO2 

4.51 

3.  19 

T1O2 

4 . 594 

2.958 

TaO  2 

4.709 

3.065 

Sn0  2 

4.737 

3.  186 

Nb02 

4.77 

2.96 

Te02 

4.79 

3.77 

WO  2 

4.86 

2.77 

M002 

4.86 

2.79 

PD02 

4.946 

3.  3'’9 

n 

1.6525 


t 

t 


2.3-2 


Ti  tends  to  form  a  system  of  nonstoch iometric  oxides 
Metallic 


TABLE  9 


RMO4  ARRANGEMENT 


Compound 

a 

c 

AlSb04 

4.510 

2.961 

CrSb04 

4.577 

3.042 

GaSb04 

4.59 

3.03 

RhSb04 

4.601 

3.100 

RhV04 

4.607 

2.923 

FeSb04 

4.623 

3.011 

CrTa04 

4.626 

3 . 009 

CrNb04 

4.635 

3.005 

FeTa04 

4.672 

3.042 

FeNb04 

4.68 

3.05 

RhTa04 

4.684 

3.020 

RbNb04 

4.686 

3.014 
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The  rutile  materials  tend  to  exhibit  high  dispersion. 

Studies  of  titanium  dioxide  films  have  been  made  by  Cherepanova 
and  Titova, 25, 26  schiller  et  al.,^"^  Joseph  and  Gagnaire^®  and 
Perney  et  al.2'^  The  index  of  refraction  of  Ti02  vary  by  10% 

depending  on  the  preparation  procedure.  The  absorption  coef¬ 
ficient,  K,  can  be  reduced  to  10~4  vvith  proper  care.  It  can  also 
be  raised  very  high  with  the  addition  of  impurities  such  as 
copper. Gagnaire  and  Joseph  also  studied  the  Ti/Ti02 
system. Samara  and  Peercy  measured  the  pressure  and  temperature 
dependence  of  c g  for  rutile. 32  optical  properties  of  the 

metals  RUO2  and  Ir02  were  studied  by  Goel,  et  al.33  while  the 
Raman  spectrum  of  Te02  was  taken  by  Pine  and  Dresse  lhaus .  34  -piig 
optical  properties  of  Sn02  are  discussed  in  references  35-37. 

In  addition  to  the  Rutile  structure,  M0O2  and  WO2  can  also 
have  a  monoclinic  structure  closely  resembling  the  Rutile  struc¬ 
ture.  Since  M0O2  and  WO2  can  occur  in  more  than  one  phase,  the 
attract  iveness  of  the  '^Ox''^l-x‘^2  system  is  considerably  reduced. 
This  consideration  is  accentuated  by  the  implications  of  a 
monoclinic  structure  for  the  anisotropy  of  the  dielectric 
constant.  The  dioxides  of  Re  and  Tc  also  have  this  monoclinic 
structure.  The  lattice  parameters  for  these  materials  are  given 
in  Table  10,  along  with  those  for  (monoclinic)  VO2,  which  has  a 
closely  related  structure.  Since  V  easily  forms  a  series  of 
nonstoichiometric  oxides  and  VO2  is  not  only  imperfectly  matched 
to  the  crystal  structure  of  other  oxides,  but  undergoes  a  metal- 
to-insulator  phase  transition  at  68  C,  it  is  not  a  promising  com¬ 
pound.  We  conclude  that  although  the  dioxides  of  a  Re,  W,  and  Mo 
are  closely  matched,  the  materials  listed  in  Table  10  are 
not  good  candidates  for  continuously  modulated  devices.  Studies 
on  tungston  oxide  films  are  reported  in  References  38-40. 

The  optical  properties  of  VO2  were  measured  by  Mokerov  and 
Begishev  as  a  function  of  temperature . 'I  ^  Vanadium  oxide  films 
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TABLE  10 


MONOCLINIC  DIOXIDES 


Compound 

a 

b 

c 

VO  2 

5.743 

4.517 

5.375 

Tc02 

R©0  T 

5.53 

4.79 

5.53 

5.562 

4.838 

5.561 

WO  2 

5.565 

4.892 

5.650 

MoO-) 

5.584 

4.842 

5.608 

were  studied  by  Gur'yanov  and  Terukov.^^  other  vanadium  oxides 
were  studied  by  Mokerov  et  al.,^3  Mokerov  and  Ryabinin , and 
Van  Hove  et.al.'^^ 

Many  dioxides  have  the  cubic  fluorite  structure;  these  are 
listed  in  Table  11.  From  Table  11  we  see  that  Cmj^Am  f 

PUj^Ce/^_j^02  f  and  U02Pt02»  should  form  continuous  solutions,  with 
lattice  matchings  of  better  than  0.3%.  Taguchi  et  al.,46 
studied  ZiO-.  films,  and  found  that  the  dielectric  constant  varied 
from  1.72  to  2.0,  depending  on  the  gas  pressure  at  which  the 
films  were  formed. 

Several  other  oxides  were  investigated  in  our  literature 
search,  none  of  which  appears  promising  for  device  applications. 
Tanaka47  gives  curves  for  several  low-index,  high-absorbance 
oxides,  which  are  listed  in  Table  12.  The  best  of  these  is  MnO, 
which  has  the  rocksalt  structure.  The  high  absorbance  of  one  of 
these  compounds,  Re^O^,  was  also  studied  by  Schlegal  et  al,48,49 
High  absorbance  was  also  studied  for  Co  impurities  in  AI2O3  by 
Niklasson  and  Grangvist.^^ 

The  tetragonol  materials,  DyV04  and  DyAsO^  are  questionable 
because  their  refractive  indices  are  low  and  very  similar. The 
indices  for  DyVO^  are  2.02  and  2.25;  for  DyAsO^  they  are  1.84  and 
1.96.  Feldman  et  al.  reported  on  the  Faraday  rotation  exhibited 
oy  Bi|^2Ge02Q,  ^^12"^^^20^  ZnO,  -2^  r  and  Ga20^»72 

These  materials  would  be  useful  in  certain  cases  due  to  their 
magneto-optic  properties.  The  dielectric  constants  of  Bij^2SiO20 
were  investigated  by  Reza  et  al.^^  Thin  films  of  Bi203  contain  a 
high  density  of  voids  but  the  refractive  index  of  about  2.55  is 
not  sensitive  to  variations  in  stoichiometry . 54 

The  remainder  of  our  literature  search  pertained  to 
materials  commonly  used  in  the  glassy  state.  Hing  investigated 
losses  in  sintered  alumina.  55  found  that  both  Mg  and  Ca , 


TABLE  11 


CUBIC  FLUORITE-STRUCTURE  OXIDES 


Compound 

a 

Zr02 

5.07 

Hf02 

5.115 

Tb02 

5.220 

Cm02 

5.372 

Ain02 

5.376 

Pu02 

5.3960 

Ce02 

5.4110 

UO2 

5.4682 

Pr02 

5.4694 

Pa02 

5.505 

ThO, 

5.5997 

which  are  often  used  to  control  sintering,  segregated  to  the 
grain  boundaries.  High  absorption  occurs  in  sintered  alumina  due 
to  residual  porosity,  grain-boundary  contamination,  and  second- 
phase  scattering.  Best  results  were  obtained  using  MgF  as  a 
dopant.  The  effect  of  stoichiometry  on  a  common  companion  of 
alumina  in  glasses,  Si02,  was  reported  by  Zuther,  et  al.  who 
studied  the  glass  system  Ga203-Ge02-As  7,05 . 56  ^he  refractive 
indices  of  the  monoclinic  materials  PbGeO^  and  PbSi03  were 
measured  by  Sugi  et  al,57  3^^  3,-q  given  in  Table  13.  These  ani¬ 
sotropic  materials  have  similar  refractive  indices  with  high 
dispersion.  The  refractive  indices  of  oxide  films  on  Ge  and  GaAs 
are  reported  respectively  by  Nazarenko  and  Rastrenenko58  and  Umeno 
et  al.58  Finally,  Ovcharenko  and  Yakhkind  looked  at  several  oxide 
additives  in  tellurite  glasses. 50  Their  results  are  presented  in 
Table  14,  which  gives  an  extrapolation  of  the  effective  indices 
of  these  additives  to  a  100%  concentration.  The  refractive  index 
of  Te02-Ba0  is  2.0976,  and  that  of  Te02~W03  is  2.1939.  We  take 
the  coexistence  of  oxides  and  tellurides  as  a  transition  to  our 
next  paragraph,  the  cha Icogen ides . 

7 . 4  Cha Icogen  i  des 

The  II /.VI  compounds  are  a  well  understood  and  well  developed 
class  of  materials  in  which  we  make  the  transition  from  oxides  to 
cha Icogen ides .  This  is  because  ZnO  and  BeO  have  the  wurtzite 
structure,  although  their  disparate  respective  lattice  constants 
of  3.25  and  2.69  angstroms  make  them  unattractive  candidates  for 
a  solute  system,  and  because  the  rocksalt  structure  contains  both 
oxides  and  cha Icogen ides .  The  wurtzite  compounds  are  presented 
in  Table  15,  and  the  rocksalt  structured  cha Icogen ides  are  given 
in  Table  16.  Copper  halides  are  also  listed  in  Table  15.  None 
of  the  wurtzite  systems  appears  attractive  for  low-defect  solid 
solutions,  due  to  the  incompatibility  of  their  lattice  parameters. 
The  same  comment  applies  to  the  rocksalt  structure  chalcogenides. 


TABLE  13 


PbGe03  and  PbSiO^ 


PbGeO- 


11.57 


7 . 32 


12.62 


2.0337 


2.0411 


2.0506 


PbSiO: 


1.947 


1.961 


1.968 


*Taken  at  579.1  nm.  These  materials  have  high  dispersion, 


TABLE  15 


WURTZITE  STRUCTURE 


Compound 

a 

c 

ea> 

(iJ 

ides 

CuH 

2. 8B3 

4.614 

CuCl 

3.91 

6.42 

C'jBr 

4.06 

6.66 

4 . 4 

Cur 

4.31 

7.09 

Agl 

4.580 

7.494 

des 

BeO 

2 .698 

4. 380 

3.0 

381 

574 

ZnO 

3.250 

5.207 

3.7 

407 

583 

Icogenides 

ZnS 

3.811 

6.234 

MnS 

3.976 

6.432 

ZnSe 

3.98 

6.53 

5.9 

MnTe 

4.087 

6.701 

MnSe 

4.  12 

6.72 

CdS 

4.  135 

6.749 

5.6 

228 

30  5 

235 

ZnTe 

4.27 

6.99 

CdSe 

4  .  30 

7.02 

6.2 

171 

214 

MgTe 

4.52 

7.  33 

98 


TABLE  16 


ROCKSALT  STRUCTURE  CHALCOGEN I DES 


Compound 

a 

MnS 

5.22 

CdS 

5.516 

5.30 

10 

CaS 

5 . 689 

6.67 

PbS 

5.9436 

17.2 

20  5 

SuSe 

6.020 

22 . 9 

PbSe 

6.  124 

280 

SnTe 


6.313 


45.0 


Another  reason  tor  passing  over  these  materials  is  that  ZnS, 
Zns-Se,  CdS,  MnS,  and  MnSe  can  assume  both  the  wurtzite  and 
zincblende  structures.  CdS  and  MnS  can  have  the  rocksalt  struc¬ 
ture  as  well.  Furthermore,  many  II/VI  materials  do  not  resist 
moisture  well. 

■Among  the  zincblende  structure  chalcogenides ,  Hg^Cdj^_j^Te 
appears  to  be  the  most  likely  candidate  system,  as  can  be  seen 
from  Table  17.  This  system  is  especially  attractive  due  to  the 
high  and  very  different  dielectric  constants  for  HgTe  and  CdTe 
respectively  of  14.0  and  7.3.  In  fact,  Hgj^Cd^_j^Te  is  known  to 
form  a  continuous  solution  as  a  function  of  x.  The  x  value  must 
be  kept  above  0.15  to  maintain  a  bandgap.  The  reflection 
spectrum  of  Hg  ^ -^^Cd  ^  2  5"^®  '^^s  taken  by  Kavalyauskas  et  al.^1 
These  authors  obtained  a  dielectric  constant  of  12  for  this  solu¬ 
tion.  The  refractive  index  of  ZnS  as  a  function  of  temperature 
is  given  by  Thompson  et  al.^2  solutions  of  cubic  CdSe^.^f 
Cdi_x^'^x^^'  and  ZnSj^Se]^_jj  were  investigated  respectively  by 
Kainthla  et  al.^2  Antoszewski  and  Pecold,^'^  and  by  Mach 
et  il.^^  Klassen  and  Ossipyan  looked  at  the  effect  of  disloca¬ 
tions  on  CdS  and  CdSe.  They  found  that  the  inhomogeneities 
oausei  by  dislocations  caused  focusing  and  channeling  of  light 
passing  through  these  materials. 


Many  transition  me ta 1/d i cha Icogen ide  compounds  take  the 
pyrite  structure  (see  Table  18).  From  Table  18  we  see  that 
7u^OS|_j^  Svolutions  with  either  S,  Se ,  or  Te  should  make  excellent 
lattice  matches.  Infrared  measurements  on  several  pyrites  were 
reported  by  Anastassakis  and  Perry.  67  verble  and  Humphrey*’® 
measured  the  infrared  and  Raman  spectra  of  MnS2^  and  Ushioda 
reported  on  Raman  measurements  on  FeS2.6^ 


One  reason  why  the  chalcogenides  can  be  interesting  is  that 
sulfur  and  selenium  have  similar  crystal  radii,  of  1.84  and  1.88 
angstroms  respectively,  so  that  several  solutions  of  Sj^Sei^^ 


are 


TABLE  17 


ZINCBLENDE  OR  SPHALERITE  STRUCTURE  CHALCOGEN IDES 


Compound 

a 

CO  -p 

BeS 

4.85 

BeSe 

5.07 

2nS 

5.409 

5.14 

8.0 

282 

BeTe 

5.54 

MnS 

5.600 

ZnSe 

5.668 

5.90 

8. 3-9.1 

207 

MnSe 

5.82 

CdS 

5.818 

5.63 

HgS 

5.85 

30.6 

HgSe 

6.084 

25.6 

ZnTe 

6.089 

7.28 

10 

178 

CdSe 

6.  16 

6.38 

HgTe 

6.429 

14.0 

CdTe 

6.480 

7.3 

140 

CuCl 

5.406 

3.6 

7.3 

155 

CuBr 

5.690 

4.4 

Cul 

6.043 

5.2 

TABLE  18 


PYRITE 


STRUCTURE 


CHALCOGENIDES 


Compound 

a 

FeS2 

5.407 

C0S2 

5.52 

RhS2 

5.574 

RuS  2 

5.59 

OSS2 

5.608 

CoSe2 

5.859 

RuSe2 

5.921 

OsSe2 

5.933 

NiSe2 

5.960 

RhSe  2 

5.99 

MnS2 

6.101 

RuTe2 

6.360 

OsTe2 

6.369 

IrTe2+ 

6.411 

RhTe2 

6.441 

21.3 


6.07-6  .  14 
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known.  In  addition,  some  chalcoyenides  have  large  dielectric 
constants.  However,  chalcogenides  tend  to  have  narrow  bandgaps, 
and  are  best  suited  for  use  in  the  infrared.  Like  the  halides, 
many  chalcogenides  form  layered  compounds  with  transition 
metals. These  materials  are  discussed  at  length  in  Wieting  and 
Schluter.^1  The  materials  are  most  commonly  grown  by  halogen 
vapor  transport .  Often  it  is  possible  to  grow  large  single 
crystals,  one  exception  being  WS2*  Many  of  these  layered  com¬ 
pounds  are  metals,  but  several  are  semiconductors.  One  of  the 
better  known  layered  structures  occurs  for  GaS  and  GaSe.  Data  on 
these  compounds  are  given  in  Table  19,  where  N  refers  to  the 
number  of  layers  within  i  inii  cell.  We  see  that  GaS  and  GaSe  do 
not  have  either  commensurate  lattice  parameters  or  particularly 
interesting  refractive  indices.  GaSe  is  an  electro-optic 
material;  its  electro-optic  properties  have  been  measured  by 
Sokolov  and  Subashiev.^3  temperature  dependence  of  the 

refractive  index  of  GaSe  was  investigated  by  Antonioli,  et  al.^4 
The  dispersion  in  GaSe  is  very  high,  with  a  refractive  index  that 
varies  from  about  2.73  at  0.5  eV  to  3.10  at  2.5  eV.74 
Reflectivity  measurements  on  GaS^Se^.^  crystals  are  presented  in 
Reference  75.  Other  optical  measurements  on  GaSe  are  given  in 
References  76-80. 

The  most  common  layered  structure  formed  by  the  transition 
metal  dichalcogenides  is  the  Cdl2  structure  shown  in  Figure  7.1a. 
Many  of  these  compounds  are  metals,  hut  except  for  SnSe2,  the 
materials  of  Table  20  are  semiconductors.  We  see  from  Table  20 
that  solutions  of  Hf/Sn/Zr  in  either  sulfur  or  selenium  are 
possibilities  for  high-index  devices  with  a  potential  for  modula¬ 
tion  of  over  30%.  The  sulphides  are  particularly  attractive, 
since  they  have  wider  bandgaps,  a  much  closer  (1%)  match  in  lat¬ 
tice  constants,  and  do  not  contain  any  metals.  Infrared  reflec¬ 
tance  measurement  on  many  of  these  compounds  were  made  by 
Lucovsky  et  al.^1  Optical  and  U.V.  measurements  on  Sn3e2  ^nd 


TABLE  19 


Gas  STRUCTURE 


Compound 

a 

c 

N 

n 

G  a  S 

3. 58S 

15.50 

4 

3.0 

GaSe 

3.746 

23.910 

6 

2.7 

I  nSe 

4. OS 

16.93 

4 

TABLE  20 


Cdr2  STRUCTURE  -  CHALCOGEN I  DBS 


Compound 

a 

c 

HfS2 

3.635 

5.837 

SnSo 

3.639 

5.868 

ZrS2 

3.662 

5.813 

SnSSe 

3.716 

6.050 

HfSe2 
ZrSe2 
SnSe  1 


3.748 

3.771 

3.811 


6.159 

6. 138-.  149 
6.137 


SnS2  were  made  by  Bertrand  et  al.^2  other  papers  on  the 
optical  properties  of  SnS2  and  SnSe2  are  given  in  references  83 
and  84.  Isomaki  and  Boehm  have  made  calculations  on  the 
dielectric  properties  of  ZrS2  • 

The  best  known  transition  metal  dichalcogenide ,  MoS2» 
belongs  to  the  octehedrally  coordinated  family  listed  in  Table 
21,  and  depicted  in  Figure  7.1b.  These  materials  have  large 
dielectric  constants  (7-24),  albeit  narrow  bandgaps.  As  in  the 
case  of  the  1T-Cdl2  materials,  we  find  that  a  match  of  lattice 
constants  separates  the  MOj^W^_^  solutions  into  sulphides  and 
selenides.  Of  these,  the  sulphides  are  harder  to  grow. 

Beal  and  Liang  have  made  reflectivity  measurements  on  WSe2 
and  WS2.^^  The  optical  properties  of  MoSe2»  which  is  relatively 
easy  to  grow,  were  measured  by  Evans  and  Hazelwood. other 
relevant  papers  are  references  88  and  89.  We  note  that  many 
transition  metal  trichalcogenides  also  take  layered  structures. 
Zwick,  et  al.  made  Raman  measurements  on  ZrS3,  ZrSe^,  ZrTe^,  and 
HfSe^.'^O  The  infrared  spectrum  of  HfSg  was  measured  by  Jandl 
and  Deslandes .  1 

The  system  Bi 2Te also  has  a  layered  structure.  It  has 
a  very  high  dielectric  constant  (20-40)  in  the  region  between  2 
and  10  urn,  where  it  has  an  extinction  coefficient  of  less  than 
one  urn.  *^2  orpiment,  HS2B3,  and  its  companion,  AS2Se3,  may 
crystallize  in  a  layered  form,  or  can  form  glasses.  Most  work  on 
these  materials  has  been  done  on  the  glass  phase.  These  glasses 
have  good  mechanical  strength  and  chemical  stability.  However, 
their  range  of  refractive  indices  is  small  (n  varies  from  2.4  to 
2.7),  and  their  ability  to  form  both  a  glass  and  a  crystalline 
phase  suggests  that  scattering  from  defects  may  be  a  problem. 
Butterf ield93  started  the  optical  properties  of  40  thin  films 
of  these  glasses.  He  attributed  the  dominant  scattering  mechan¬ 
ism  to  void  and  oxygen  impurities.  Kanchiev  and  Kokorina  concur 


TABLE  21 


MOLYBDENITE  STRUCTURE 


Compound 

a 

c 

^  OD 

M0S9 

3.16 

12.30 

10  (2.2um) 

WS2 

3.15 

12.36 

7-8 

MoSe2 

3.29 

12.93 

24 

WSeo 

3.29 

12-97 

7.3 

that  oxygen  impurities  are  a  serious  problem,  as  is  water. 94 
They  suggest  distillation  as  a  means  of  removing  these  impuri¬ 
ties.  Young95,96  finds  that  the  absorption  and  reflectivity  of 
AS2S3  films  depend  on  surface  preparation,  especially  the  finish 
and  the  temperature  at  which  the  films  were  deposited.  Other 
papers  on  the  optical  properties  of  these  glasses  are  given  in 
references  97-107.  GeS2  can  also  take  both  a  glassy  and  a 
crystalline  form,  and  is  sometimes  used  to  form  glasses  and 
AS2S3.  Raman  studies  on  GeS2  have  been  made  by  Nemanich. 107 , 108 
The  single  sulphides  and  selenides  of  Ge  and  Sn  have  an 
orthorhombic  structure  (Table  22).  These  materials  do  not  appear 
to  be  well  suited  for  continuously  modulated  optical  devices. 
Their  optical  properties  are  reported  in  references  109-116. 
Miscellaneous  references  on  chalcogenide  compounds  are  given  in 
references  117-121. 

Our  remaining  structure  for  chalcogenide  compounds  is  the 
chalcopyrite  structure.  This  structure  is  shared  by  some  arsenic 
compounds  and  thus  forms  a  bridge  with  the  Group  V  materials. 

The  chalcopyrites ,  listed  in  Table  23,  are  birefringent  and 
difficult  to  grow  in  large  crystals.  From  Table  23  we  see  that 
there  are  four  closely  matched  systems,  all  of  which  contain 
either  Cu  or  Ag ,  and  most  of  which  contain  either  A1  or  Ga.  Most 
of  these  materials  are  selenides.  Of  the  four  systems, 

CuAl^Ga 3_xS2  f  CuAl^Ga3_jjSe2  f  have  very  closely  matched  lattice 
parameters.  The  In  (  CuSe2  )  { AgS2  )  ^nd  Al  ( AgSe2  )  ^  (  CuTe2  ) 
systems  are  likely  to  form  localized  compounds  such  as  InCuS2 

or  AlCuSe2»  which  would  have  dissimilar  lattice  parameters.  The 
Raman  and  infrared  spectra  of  AgGaS2  was  taken  by  Holah  et 
al.l22  AgGaS2  has  a  direct  band  gap  of  2.75  eV,  and  a  large 
nonlinear  optical  coe f f ic lent . ^ 2 3  Many  chalcopyrites  also  take 
the  zinc-blend  structure.  The  (CuInSe2)i_x  ~  (ZZnSe)j^  system  was 
studied  by  Gan  et  al.,  in  this  context.  124  reflectivity  of 

a  number  of  chalcopyrites  was  measured  by  Rife  et  al.l25 


TABLE  22 
GeSe  STRUCTURE 


Compound 

a 

b 

c 

GeSe 

4.40 

11.82 

3.85 

SnS 

4.33 

11.18 

3.98 

GeS 

4.30 

10.44 

3.65 

SnSe 

4.46 

11.57 

4.19 

PbSnS2 

4.23 

11.42 

4.09 

TABLE  23 


COMPOUND 

C  u  A  L  S  2 
CuGaS  2 

Cu I nS  2 

CuT1S2 

CuAlSe2 
CuGaSe  2 

ZnS  iAS2 

AgFeS2 

ZnGeAs  2 

AgAlS2 

AgGaS 

Cu I nSe  2 
Ag  rnS2 
CuT  1  Sa  2 

CdGaAs  2 

AgAlSe2 

CuAlTa2 

AgGaSe2 

C  uGaT  6  2 

Ag  InSe2 

CdSnAs2 

Cu I nTe2 

AgG  aTa  2 

Ag  InTe2 


CHALCOPYRITES 


a 

5.312 

5.349 

5.517 

5.580 

5.606 

5.607 

5.608 
5.66 
5.670 
5.695 
5.743 

5.773 

5.816 

5.832 

9.942 

5.956 

5.964 

5.973 

5.994 
6.090 
6.092 
6  .  167 
6 . 288 


c 

10.42 

10.47 

11.06 

11.17 

10 . 90 
10.99 

10 . 89 

10.30 

11.153 

10 . 26 

10.26 

11.55 

11.17 

11.63 

11.224 

10.75 

11.78 

10.88 

11.91 
11.67 

1 1 .922 

12.34 

11.94 


6.406 


12.56 


7.5  Group  V  Compounds 


The  group  V  materials  have  strongly  covalent  bonds,  and  in 
general  exhibit  narrow  bandgaps.  Consequently  they  are  of 
interest  primarily  for  infrared  devices.  The  nitrogen  compounds 
tend  to  have  bond  lengths  that  are  not  well  matched  to  those  of 
P,  As,  or  Sb,  so  we  will  concentrate  our  attention  on  the  latter 
chemicals,  and  defer  a  discussion  of  nitrogen  compounds  to  the 
end  of  this  paragraph. 


By  far  the  most  well  developed  and  most  commonly  used  group 
V  materials  are  the  III-V  compounds  with  the  zincblende 
structure.  Simply  by  virtue  of  current  materials  processing 
capability,  these  compounds  deserve  attention,  and  are  listed  in 
Table  24,  in  which  are  included  the  closely  related  materials  Si 
and  Ge.  These  materials  do  suffer  from  generally  narrow 
bandgaps  and  low  melting  temperatures.  However,  they  can  be 
grown  on  a  large  scale  with  a  very  low  level-  of  defects.  From 
Table  24  we  see  that  an  Al^Gaj^.j^  solution  with  either  P,  As,  or 
Sb  will  give  a  good  lattice  match,  albeit  with  a  variation  in 
dielectic  constant  of  only  about  20%.  Continuous  modulation  of 
these  materials  should  reduce  the  bandbending  associated  with 
interfaces .* 26  thus  seems  that  these  materials  are  attractive 

for  narrow-band  ignition  filters,  provided  that  input  power 
levels  are  not  too  high.  References  127-130  concern  the  optical 
properties  of  some  of  these  materials. 

Our  final  group  of  materials  are  the  nitrides  with  the  ZnO 
structure  (see  Table  25).  We  expect  these  nitrides  to  have  wider 
bandagaps  and  higher  melting  points  than  the  zincblende  compounds 
of  Table  24.  It  appears  that  the  Ga^Alj.^N  and  Nbj^Taj^.^N  systems 
are  good  candidates  for  optical  devices. 


.*■  V  V  V  ] 


TABLE  24 


III/V  COMPOUNDS  -  ZINCBLENDE  STRUCTURE 


IMPOUND 


iC 

D 


aP 

IP 


lAs 

aAs 


tP 

tAs 

aSb 

ISb 

-!.Sb 


a 

.  6 

n 

5 

00 

4.348 

4.538 

6.7 

10 . 2 
10.2 

5.4307 

3.5 

11.7 

5.4505 

5.451 

3.46 

2.94 

8.5 

7.6 

10.18 

5.62 

5.6537 

5.6574 

4.016 

9.0 

10.9 

16 

13.13 

15.8 

5.8687 

3.44 

3.08 

9.6 

12.37 

6.036 

4.  17 

3.46 

12.3 

14.6 

6.118 

6.1347 

4.24 

3.824 

14.4 

10.2 

15.7 

11.2 

6.478 


15.6 


17.88 


TABLE  25 


ZnO  STRUCTURE  NITRIDES 


COMPOUND 

a 

c 

e® 

(jO 

AIN 

3.111 

4.978 

4.8 

622 

GaN 

3.  180 

5.166 

5.8 

InN 

3.533 

5.693 

NbN 

3.017 

5 . 58  0 

TaN 

3.05 

4.94 
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CONCLUSION 


We  now  present  the  conclusions  resulting  from  our  literature 
survey  and  preliminary  analysis.  Although  he  has  not  definiti¬ 
vely  proved  all  of  these  conclusions,  the  author  believes  that 
they  are  valid  over  sufficiently  general  conditions  to  warrant 
their  entry  in  this  report.  We  assert  that; 

•  Stop  and  pass  bands  will  be  formed  whenever  the  suscep¬ 
tibility  is  periodically  modulated.'  In  fact  a  series  of  such 
bands  will  occur  as  a  function  of  the  wave  frequency.  The  band 
structure  will  become  more  sharply  defined  as  the  number  of 
repeat  regions,  or  unit  cells,  increases,  and  as  the  medium 
approaches  absolute  periodicity.  It  is  likely  that  the  sharpness 
of  the  band  edges  will  be  enhanced  if  the  imaginary  part  of  the 
material  susceptance  is  minimized. 

•  The  centers  of  successive  stop  bands  occur  at  wavelengths 
for  which  the  optical  depth  of  the  medium  taken  over  a  unit  cell 
distance  is  an  integral  number  of  half  wavelengths,  i.e.,  the 
phase  angle  is  equal  to  tt  ,  or 

d 

k/dxn(x)=Tr  .  (8.1) 

o 

o 

Thus  the  location  of  the  gaps,  or  stop  bands  is  determined  pri¬ 
marily  by  the  dimension  of  the  unit  cell,  d,  and  the  mean  value 

of  the  refractive  index,  n(x),  rather  than  the  shape  of  the  modu¬ 
lation  profile. 

•  The  width  of  the  band  gaps,  or  stop  bands  increases  as  the 
relative  modulation  in  the  refractive  index,  An/ri,  increases. 

For  small  modulations,  in  most  cases  the  widths  of  stop  bands  of 
a  given  order  are  roughly  proportional  to  components  of  the 
Fourier  transform  of  the  modulation  profile  of  similar  order. 


•  The  viegree  of  attenuation  of  a  wave  with  a  frequency  within 
a  stop  band  increases  toward  the  middle  of  the  band.  The  airvaunt 
of  attenuation  is  proportional  to  An/n.  If  the  number  of  repeat 
units  (unit  cells)  in  the  filter  is  sufficiently  large,  such  a 
wave  will  not  propagate  through  the  filter  at  all,  it  will  be 
totally  stopped.  The  attenuation  of  the  wave  is  not  due  to 
absorption,  it  is  due  to  reflection,  i.e.,  the  stop  bands 
correspond  to  reflection  bands. 

•  Complex  modulation  profiles  in  the  medium  susceptance  will 
lead  to  "band  spl itt ing"--the  formation  of  new  stop  bands  within 
what  would  otherwise  be  pass  bands.  Consequently  several  stop 
bands  can  be  located  at  arbitrary  frequencies  within  say,  one 
octave  of  the  first  stop  band  by  appropriate  design  of  the  modu¬ 
lation  profile.  The  locations  and  widths  of  the  various  stop 
bands  will  be  related  to  the  filter  parameters  in  the  same  manner 
as  for  a  single  stop  band.  If  it  is  desired  to  locate  stop  bands 
at  arbitrary  intervals  over  a  multioctave  region,  then  it  will 
probably  be  necessary  to  put  up  with  additional  "ghost"  stop  bands. 

•  As  7i  increases,  the  acceptance  angle  for  the  filter  will 
also  increase. 

•  Unless  performance  is  desired  over  many  octaves,  there  is  no 
apparent  reason  to  favor  one  type  of  unit  length  modulation  pro¬ 
file  (e.g.,  a  step  function  as  a  sinusoidal,  parabolic,  or 
Gaussian  profile)  over  another.  Details  of  the  modulation  pro¬ 
file  only  become  significant  in  the  upper  level  stop,  or  reflec¬ 
tion  bands.  There  may  be  reasons  based  on  materials  for  choosing 
a  given  type  of  modulation  profile.  For  example,  it  may  be  that 
impurities  will  aggregate  at  discontinuities.  However,  our  ana¬ 
lysis  in  Paragraph  2.2  suggests  that  even  if  such  an  aggregation 
were  to  occur,  it  would  not  affect  the  amount  of  absorption  in  most 
cases.  Some  effect  might  be  discernable  at  stop  bands. 


In  addition,  we  have  made  a  brief  preliminary  survey  of 
possible  materials  for  graded-index  optical  devices.  In  our 


analysis  we  assumed  that  a  minimization  of  absorption  was  para¬ 
mount,  and  accordingly  sought  systems  that  offered  a  hope  of 
achieving  a  low  density  of  defects  and  an  infrared  spectrum  rela¬ 
tively  free  of  absorption  bands.  This  implied  a  choice  of  opti¬ 
cal  materials  with  simple  molecular  formulas,  materials  that 
crystallize  with  identical  structures  and  closely  matched  lattice 
parameters,  and  materials  that  are  known  to  crystallize  in  only 
one  form.  In  particular,  we  noted  that  many  potential  candidate 
materials  are  layered  compounds.  It  may  be  possible  to  use 
layered  materials  with  dissimilar  chemical  formulas  in  adjacent 
layers  without  introducing  serious  problems  with  interface 
states . 

8 .  1  Comments 

The  author  is  confident  that  an  analytic  description  of  a 
filter  behavior  similar  to  that  discussed  in  Sections  3-6  can  be 
extended  to  complex  susceptibilities  and  three-dimensional 
systems.  Similar  problems  have  been  treated  in  three  dimensions 
in  band  theory,  lattice  dynamics,  and  the  description  of 
microwave  components.  Step  function  modulations  are  the  easiest 
to  analyze  because  the  solutions  are  in  terms  of  sines  and 
cosines,  which  are  familiar,  and  can  often  be  manipulated  into 
compact  expressions.  The  analysis  of  a  Rugate  filter  is  much 
less  tractable.  However,  in  many  cases  of  interest,  such  as 
sinusoidal  modulation,  closed  form  solutions  should  be 
obtainable.  Indeed,  solutions  to  the  Mathieu  problem  have  been 
known  for  over  half  a  century. 

Commonly,  continuously  modulated  filters  are  approximated  in 
manufacture  by  stacking  together  a  large  number  of  very  thin 
layers,  with  flat  index  profiles  across  each  layer  and  sharp 
discontinuities  between  layers.  Our  analysis  suggests  that  the 
merit  of  this  procedure  is  at  best  moot.  This  practice  will 
indeed  ]ive  a  reasonable  approximation  of  an  arbitrary  index  pro- 


file  over  3  liniteo  (one  or  two  octave)  frequency  range. 

However,  it  introduces  a  very  large  number  of  discontinuities. 

If  these  discontinuities  do,  in  fact,  result  in  deterioration  of 
the  material  properties,  then  these  discontinuities  are  self 
defeating.  They  will  also  affect  the  structure  of  high-order 
bands.  Therefore  if  the  higher-order  band  structure  is  impor¬ 
tant,  the  use  of  thin  layers  will  interfere  with  these  upper- 
level  bands. 

We  note  that  both  the  width  and  the  degree  of  attenuation 
or  reflection  in  the  stop  or  reflection  bands  scales  as  in,/?i. 

This  creates  difficulties  if  notch  filters  are  desired.  It 
appears  that  the  primary  way  out  of  this  dilemma  is  to  bu i 1 i  a 
filter  with  many  repeat  distances,  i.e.,  a  very  long  filter.  It 
may  also  be  possible  to  achieve  some  design  freedom  by 
prescribing  profiles  that  have  narrow  potential  wells  (regions  of 
large  susceptibility)  separated  by  relatively  wide  distances, 

8 . 2  Implications  and  Suggestions 

On  the  basis  of  this  interim  survey  the  author  concludes 
that  the  choice  of  optimum  filter  materials  and  manufacturing 
should  be  ciominated  by  materials  and  manu  f  ac  tu  r  i  ng  questions 
ratner  than  by  filter  design  criteria.  The  primary  concern 
should  be  to  minimize  intrinsic  losses.  This  can  be  achieved  by 
selecting  materials  whose  imaginary  component  of  the  refractive 
index,  nj^,  is  small  compared  to  the  real  component,  nj_,  and  by 
minimizing  local  inhomogeneities  which  will  add  to  losses,  as 
described  in  Paragraph  2.2.  The  second  criteria  should  be  to 
select  a  material  that  is  immune  to  intense  levels  of  signal 
energy.  This  implies  a  low  variation  of  n(x)  with  temperature, 
and  a  high  melting  or  phase  transition  temperature.  Comments  on 
nonlinear  susceptibilities  are  beyond  the  scope  of  this  report. 


In  general,  it  will  be  desirable  to  adjust  the  widths  of  the 
stop  or  reflection  bands  and  to  create  multiple  stop  bands  at 


arbitrary  f  r-^quenc  i  es .  This  will  require  tailor-nade  refraction 
indices,  which  are  achieved  by  appropriate  solid  solutions.  If 
i  nhomoqene i t ies  in  the  resulting  material  are  to  be  minimized, 
the  constituents  o[  the  solution  should  be  as  compatible 
(mutually  soluble)  as  possible.  This  suggests  similar  struc¬ 
tures,  chemical  compositions,  and  bond  lengths. 

In  summary,  we  believe  that  optimum  materials  should  have 
the  following  properties: 

•  nj^  <<  n^. 

•  Amenability  to  manufacture  with  a  low  concentration  of 
defects 

•  dn/dT  -*•  0 

•  A  high  phase  transition  temperature 

•  Good  solubility  with  other  materials  possessing  a  simi¬ 
lar  structure  but  as  different  a  refractive  index  as 
poss  ible . 

This  interim  report  is  by  its  very  nature  incomplete.  If  we 
are  to  obtain  a  reasonably  complete  analytic  picture  of  maximum 
filter  design,  several  steps  should  be  undertaken.  We  need  to 
make  analytic  calculations  of  the  dependence  of  the  imaginary 
part  of  the  wavenumber,  kj^,  and  the  reflection  bandwidth,  Aw  on 
the  relative  modulation  depth,  An/rT.  These  calculations  should 
be  made  for  both  the  Kron ig-Penney  and  the  Mathieu  problems,  and  the 
respective  results  should  be  compared.  The  analytic  calculations 
should  be  augmented  by  numerical  calculations,  and  should  be 
extended  to  include  multiply  periodic  filters.  We  should  also 
make  numerical  comparisons  between  the  Kron ig-Penney  and  the  Mathieu 
contours  in  regard  to  their  higher-order  band  structure  and  the 
side  lobes  in  their  reflection  bands. 


We  shoul'i  extend  our  calculation  to  include  complex  suscep¬ 
tibilities,  including  both  analytic  and  numerical  work. 
Formulations  should  be  given  for  additional  profiles,  e.g.,  as 
sawtooth  and  parabolic  profiles.  In  particular,  we  should 
investigate  the  application  of  a  Kramers-Kron ig  analysis  to  a 
3ragg  filter.  The  Kramers-Kron ig  relation  is  based  solely  upon 
causality.  Therefore  it  is  reasonable  to  assume  that  it  will 
apply  to  the  filters  in  question.  However,  since  we  get  stop 
bands  even  with  ideal  materials  that  have  purely  real  suscep¬ 
tibilities,  the  applicability  of  a  Kramers-Kron ig  analysis  is  not 
immediately  obvious.  Assuming  that  it  does  apply,  it  would  be 
very  useful  both  in  filter  design  and  in  reducing  data  on  filters. 

Finally,  our  analytic  formulations  should  be  extended  to 
three  dimensions,  particularly  for  the  case  of  the  Mathieu  problem 
Such  calculations  can  most  readily  be  made  in  Cartesian  coor¬ 
dinates,  but  calculations  in  a  cylindrical  coordinate  system 
should  also  be  attempted.  If  time  permits,  an  analytic  treatment 
of  the  effect  of  terminations  and  interfaces  should  also  be  made. 


APPENDIX  A 


MATHEMATICAL  UNDERPINNINGS  OF  WAVE  PROPAGATION  IN  PERIODIC  MEDIA 

In  this  appendix  we  derive  some  o£  the  equations  which  form 
a  basis  for  the  discussion  of  Section  3.  The  fundamental  pro¬ 
perty  that  facilitates  the  description  of  this  problem  is  that 
the  waves  obey  a  second-order  equation  of  motion, 

■4,"  fix)  ^  ~  0  ,  (A.l) 


in  which  no 
is  that  the 
the  path  of 
be  involved 
many  fields 
vibrations, 
band  theory 


first  derivative  term  appears.  The  second  property 
response  of  the  media  to  the  wave  is  periodic  along 
the  wave.  No  other  properties  of  the  system  need  to 
in  our  discussion.  Consequently  our  remarks  apply  to 
other  than  optics,  including  acoustics,  lattice 
transmission  line  theory,  waveguide  theory,  and  the 
of  solids.  In  Section  3  we  use  . 


f(x)  =  k^n^(x)  =  k^e(x)  =  k^(x) 
o  o 

A. 1  Form  for  Transfer  Matrix  Elements 


(A. 2) 


A  second-order  differential  equation  will  have  two  indepen¬ 
dent  solutions  for  a  given  set  of  parameters  (e.g.,  k  and  w). 

Let  these  be  41  and  ^2'  Then,  from  Equation  (A.l), 


(0  =  -  f  (  X)  <{.  ^<{1  2 

and 


(A.  3) 


(A.  4) 


Subtracting  (A. 4)  from  (A. 3)  we  have 


This  is  a  key  result,  which  as  we  see,  stems  from  the  form  of 
{ A. 1 ) .  Now 


dx  ^*1^2  ~  2*^ 


4>i  -  tot; 


(  A .  6  ) 


Consequently,  we  find 


d7  ^*1*^2  ■  °  ' 


(A. 7) 


j^t2  ~  t2t[  =  constant  . 


(A. 8) 


The  quantity  in  the  l.h.s.  of  (A. 8)  is  known  as  the  Wronskian  for 
(A.l).  The  critical  point  is  that  the  Wronskian  is  independent  of 
the  argument  of  (f  and  <^2'  long  as  they  are  taken  at  the  same 
point.  Note  that  if  we  had  an  equation  of  the  form 


({/"  +  y(x)i|;'  +  f(x)  i|;  =  0  , 


we  would  obtain  the  result 


(A. 9) 


1^2  ”  ‘^2'*'1  ^  {-/Y(x)dx) 


(A. 10) 


We  would  need  to  use  (A. 10)  instead  of  (A. 8)  in  regard  to 
the  full  wave  Equations  (4.1)  -  (4.4). 

We  are  now  in  a  position  to  derive  the  expressions  for  the 
matrix  elements  Equations  (3.17)  -  (3.21).  We  begin  with  a 
general  solution  to  (A.l),  which  is  a  linear  combination  of  the 
independent  solutions  ifi  ^  and  ^2' 


Y  +  C24>  2 


(A. 11) 


Differentiation  of  (A. 11)  yields 


>)'  ’  =  c.<t>  ;  +  I 


(A. 12) 


(A. 11)  and  (A.12): 

'^l^*^l'^2  ~  '^2'^P  '  (A. 13) 

D 

^2^^1^2  ■  *  (A. 14) 

1.  -  --^ 

D 

Since  the  c  coefficients  are  the  quotients  of  Wronskians  of  (A.l), 
they  are  indeed  constants,  by  (A. 8)’.  We  may  therefore  use  the 
argument  at  x'  to  for  and  C2  in  (A. 11): 

\p(x)  =  Ci(x'  )4i  ^(x)  +C2(x’)(j)2(x)  .  (A. 15) 

Substituting  from  (A. 13)  and  (A. 14)  in  (A. 15),  we  find 

\l;(x)  =-^[ijj(x‘)<()2(x')  -i()'(x)())2(x’)]()>j^(x) 

+  ^  [^Mx' )^^(x' )  -  i}»(x' )♦  J(x' )]4.2(x)  .  (A. 16) 

We  may  rearrange  the  terms  in  (A. 16)  to  obtain  an  expression  for 
ii{x)  in  terms  of  4'(x')  and  i()'(x'); 

it<(x)  =-^[((i^(x)(ii^(x')  -4i2(x)4i^{x')]i|»(x') 

+  "5  ~  ()»j^(x)4i2(x')]<^’(x') 

=  b^^i};(x')  +  O.E.D.  (A. 17) 

A  similar  argument  will  give  the  expressions  for  the  matrix  ele- 


We  may  solve  for  c^  and  C2  in 


>)'*2  ■  't'  '^2 


-  I)*  <t> 


rnGHts  ^21  ^22* 


A . 2  Dispersion  Relation 

We  now  consider  a  different  point  of  view,  in  which  we 
express  ( x )  in  terns  of  its  value  at  x' ,  without  regard  to 
V  ' ( X '  ) .  We  again  use  ( A .  11) ,  and  wr  i  te 

t);  (  X  )  =  X  (  X  ,  X '  )ij)  (  X  '  )  (A. 18) 

Using  (A. 13),  and  letting  '1;  ( x '  )  also  be  of  the  form  {A. 11),  there 
will  be  some  matrix  A  with  elements-  a ,  such  that 

(j)j^(x)=a]^2^(<)j^(x')+a]^2<)>2^’^’^  (A. 19) 


a  nd 


;j>2(x)  -  a2]^'{i]^(x')  +  32^24*2^^*^  *  (A. 20) 

Equations  (A. 11)  and  ( A . 1 8 ) - ( A. 20 )  give  us 

i|j(x)  =  (^l^n  +  ^2^2J  4'i(x')  +  (^1312  ^2^22^  <|)2(x’) 

=  X  [  C2^<}i  ^  (  X  '  )  +  C2<ti  2  ( H  •  (A.  21) 

From  (A. 21)  we  obtain  the  matrix  eigenvalue  equation 


ail  ■  ^ 

®12 
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^22  ~ 

=  0 


(A. 22) 


or 


X  -  +  a22)^  +  ^'^11^22  "  ^12  ^21^  "  ° 


( A. 23 ) 


We  next  show  that  A  =  1.  Note  however,  that  since  A  =  1 ,  we  find 

X^  -  (a^^  -f-  a22)^  +  1  =  0  ,  (A. 24) 


or 
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Tj.[Al  =  X  +  X 


-1 


(  A .  2  5  ) 


We  now  show  that  A  =  1.  We  refer  to  the  (constant) 

Wronskian  D.  Using  (A. 19)  and  (A. 20)  we  write 

D  =  4>  j^(x)ij)^(x)  -  *  J^{x)4>  ^(x) 

-  '  (A. 26) 


or 


°  ^  ^^11^22  "  ^12^21^^'^!^’'' ^'^2^^' ^  ~  ij>|(x' )i}>2(’‘’ )] 

i .  e.  , 


There  fore 


D  =  AD 


A  =  1 


(A. 27) 


(A. 28) 


(A. 29) 


The  results  we  have  derived  have  depended  on  the  form  (A.l) 
for  the  wave  equation,  and  in  particular  on  the  consequence  (A. 8) 
However,  the  expression  (A.  18)  is  most  useful  when  ij;  is  periodic 
in  X.  This  will  be  the  case  when  f(x)  is  periodic,  as  we  now 
show.  We  may  construct  a  function  A(x)  such  that 


i|;(x)  =  A(x)  eikx 

Substituting  (A. 29)  in  (A.l)  we  find 

A"  +  2ikA'  +  [f(x)  -  k2j  a  =  0 


(A. 30 ) 


(A. 31) 


If  f(x)  is  periodic  with  period  d,  then  from  (A. 31)  A  must 
be  also.  It  is  often  useful  to  rewrite  (A. 30)  in  another  form. 
We  let 


X  -  x'  =  d  , 


(A. 32) 
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